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Ywo  approaches  are  used  to  explore  the  effects  of  shear  on  short-wavelength  internal 
waves.  In  the  first,  the  Taylor-Goldstein  equation  is  solved  exactly.-  The  solutions  show  that 
the  amplitude  of  partially  reflected  waves  is  exponentially  small  as  c[U<y,  where  c  is  the  hor¬ 
izontal  phase  speed  of  the  short  waves  and  Uo  is  the  maximum  mean  flow  velocity: v  The  solu¬ 
tions  also  reveal  the  inaccuracy  of  WKB  predictions  when  applied  to  a  curved  velocity  profile 
with  a  minimum  Richardson  number  of  order  unity. 

To  investigate  internal  waves  in  an  inertial  current,  ray  calculations  are  made.  This 
second  approach  reveals  that  the  process  of  refractive  convergence,  which  includes  the  critical- 
layer  interaction  as  a  special  case,  operates  at  virtually  all  phases  of  the  inertial  oscillation  and 
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'  affects  short  waves  of  nearly  all  frequencies,  it  is  also  found,  in  contrast  to  the  results  of  steady 
shear  analyses,  that  short  waves  with  phase  speeds  less  that  the  mean  flow  maximum  can  pro¬ 
pagate  for  several  inertial  periods  without  becoming  unstable,  and  conversely  that  waves  with 
initial  phase  speeds  of  two  or  three  times  the  mean  flow  maximum  can  quickly  become 
focussed  to  unstably  high  amplitudes. 

The  final  section  examines  the  mean  flow  induced  by  three-dimensional,  low- 
frequency,  internal  wave  packets.  Rotation  alters  the  character  of  the  flow  so  that  the  mean 
momentum  is  not  equal  to  E/c,  where.  E  is  the  intrinsic  energy  density.  The  generation  of 
inertial  waves  by  the  internal  wave  field,  as  found  by  Hasselmann  (1970)  for  a  wave  field  that  is 
statistically  homogeneous  in  the  horizontal,  is  not  predicted  by  a  calculation  that  incorporates 
horizontal  variations., 


/; 
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CHAPTER  l 


INTRODUCTION 

A  challenging  feature  of  the  internal  wave  held  is  its  maintenance  of  high  and  uniform 
energy  levels.  A  central  difficulty  arises  from  our  poor  understanding  of  the  mechanisms  by 
which  internal  waves  dissipate.  To  understand  such  mechanisms  is  to  understand  that  portion 
of  the  spectrum  where  dissipation  is  most  effective:  the  short  wavelengths. 

Several  factors  complicate  the  analysis  of  short  waves.  Because  of  their  small  phase 
velocities,  short  waves  are  strongly  refracted  by  background  shears.  In  the  deep  ocean,  these 
shears  are  produced  predominantly  by  inertial  oscillations  and  therefore  vary  in  time  as  well  as 
in  space.  Short  waves  are  also  sensitive  to  finestructure  in  the  mean  density  and  velocity 
profiles  (Mysak,  1976),  and  studies  of  weakly  nonlinear  waves  predict  such  rapid  time  scales 
for  the  interaction  of  short  waves  that  according  to  Holloway  (1980)  they  "contradict  their 
premises". 

No  single  theory  can  handle  all  of  these  processes.  In  this  thesis,  we  concentrate  on 
one:  the  effects  of  shear,  particularly  inertial  shear,  on  the  propagation  of  short  waves.  The 
objective  is  to  explore  the  mechanisms  by  which  inertial  currents  refract  short  waves  to 
unstable  amplitudes.  The  major  results  are  derived  from  ray-tracing  calculations  (Chapter  5), 
with  Chapters  3  and  4  devoted  to  assessing  the  accuracy  of  ray  theory  over  the  range  of  condi¬ 
tions  typically  found  in  the  deep  ocean. 

One  dissipation  mechanism  to  be  examined  is  the  critical-layer  interaction.  A  critical 
layer  is  defined  as  the  depth  at  which  the  horizontal,  Doppler-shifted  phase  velocity  equals  the 
mean  flow  velocity.  Within  the  WKB  approximation,  the  critical-layer  interaction  appears  as  a 
pattern  of  converging  rays  (figure  1.1).  There  is  good  reason  to  suspect  this  critical-layer  con¬ 
vergence  as  important  for  short-wave  dynamics  and  dissipation.  Many  studies  suggest  that  a 
narrow  critical-layer  region  that  remains  clearly  defined  is  connected  with  overturning  and 
mean-flow  accelerations  (Booker  &  Bretherton,  1967  (hereafter  BB);  Thorpe,  1981:  Koop, 
1981;  Fritts,  1978),  and  since  short-wave  phase  speeds  typically  fall  within  the  range  of  the 
currents  through  which  they  propagate  (Sanford,  1975;  Pinkel,  1981),  critical-layer  conditions 
are  frequently  satisfied. 

However,  for  inertial  shears,  we  show  in  Chapter  5  that  convergence  develops  for 
short  waves  with  phase  speeds  two  or  three  times  the  mean-flow  maximum  and  that  conver¬ 
gence  is  equally  likely  when  critical-layer  conditions  are  not  satisfied. 


Vt. 


Figure  1.1.  Ray  trajectories  for  short  waves  in  an  inertial  current.  The  x— component  of  the 
inertial  velocity  is  given  by  i/-C/osin/3zsin/f,  where  /3— 2ir/\,  and  /— 2w/T,  are  the  inertial 
wavenumber  and  frequency  respectively.  The  velocity  U  is  identically  zero  at  r-0  and  r-.5  T, 
and  attains  the  profile  shown  on  the  vertical  axis  at  r— .25  7^.  Behavior  resembling  the  critical- 
layer  simulations  of  Thorpe  (1981)  begins  near  r-.l  7)  and  z-l.OX,,  where  rays  merge. 

Rays  terminate  when  the  intrinsic  energy  density  E  (Chapter  5)  surpasses  50  times  its  ini¬ 
tial  value  or  when  the  Richardson  number  drops  below  one-quarter,  conditions  which  arise  in 
more  than  half  of  the  trajectories.  Initial  conditions  include  (see  Chapter  5  for  notation  and 
assumptions)  w'/Af-. 2,  c/(/o-1.0,  and  7-iV2/(/32C/o2)-4.0. 


These  results  are  unlike  those  obtained  from  the  Taylor-Goldstein  equation  used  by 
BB  and  in  many  other  critical-layer  studies.  The  Taylor-Goldstein  model  (or  its  extension  to 
rotating  fluids  (Jones,  1967) )  is  inappropriate  for  the  description  of  short  waves  in  the  deep 
ocean  for  at  least  two  reasons.  First,  it  is  based  on  the  assumption  of  steady  shears,  but  as 
figure  1.1  illustrates,  the  time  dependence  associated  with  inertial  currents  affects  the  propaga¬ 
tion  of  short  waves  within  one-tenth  of  an  inertial  period.  Secondly,  the  Taylor-Goldstein 
equation  predicts  a  distribution  of  wave  amplitude  and  intrinsic  frequency  in  the  vertical  that 
applies  to,  for  example,  waves  produced  by  steady  flow  over  topography.  Short  waves  in  the 
ocean’s  interior,  however,  may  be  aligned  differently,  depending  in  part  on  their  generation 
mechanism  and  in  part  on  the  way  they  are  refracted  by  time  dependent  currents.  An  example 
is  the  case  considered  by  Hartman  (1975)  and  extended  to  time  dependent  flows  by  Thorpe 
(1978,  Appendix).  Here,  the  spread  of  "critical  layers"  over  a  continuous  range  of  depths  leads 
to  behavior  contrasting  the  BB  results. 

Before  examining  the  effects  of  inertial  currents  on  short  waves,  a  review  of  the  BB 
and  Hartman  analyses  and  a  general  discussion  of  critical  layers  in  steady  flows  is  presented  in 
Chapter  2.  Then,  in  preparation  for  Chapter  5,  we  estimate  the  validity  of  the  slowly  varying 
approximation  necessary  for  ray  theory.  This  is  begun  in  Chapter  3  by  describing  the  scales  of 
short  waves  in  the  deep  ocean  and  by  identifying  a  criterion  for  WKB  accuracy.  In  Chapter  4 
the  magnitude  of  the  energy  reflected  from  strong  (but  steady)  shears  is  assessed  as  a  possible 
source  for  the  failure  of  the  WKB  approximation. 

In  Chapter  5,  ray  theory  is  used  to  describe  the  propagation  of  short  waves  in  inertial 
shears.  The  results  are  contrasted  with  the  steady-shear  predictions  reviewed  in  Chapter  2. 
The  effect  of  dispersion  and  geometrical  spreading  on  the  convergence  of  short-wave  packets  is 
considered  in  Chapter  5,  Section  C.  Chapter  6  concludes  with  a  summary  and  discussion  of 
results. 

Pans  of  the  derivations  in  Chapter  4  and  Chapter  5,  Section  C,  are  left  to  Appendices 
A  and  B.  In  Appendix  C,  which  contains  its  own  acknowledgements  and  references,  we  exam¬ 
ine  the  mean  flow  and  wake  produced  by  a  low-frequency  wave  packet. 

Throughout  it  is  assumed  that  the  Vaisala  frequency  is  constant  and  that  the  back¬ 
ground  current  is  horizontal  with  vertical  shear,  as  is  produced  by  an  inertial  oscillation.  Time 
dependence  and  rotation  are  neglected  until  Chapter  5.  The  terms  background  flow  and  mean 
flow  are  used  interchangeably  even  though  the  mean  flow  may  itself  be  an  internal  wave.  The 
mean  flow  also  may  have  a  frequency  and  vertical  wavenumber  not  much  less  than  the  short 
waves  it  is  refracting  but  will  always  be  distinguishable  by  its  larger  horizontal  particle  velocities 
and  infinite  horizontal  scale. 
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The  major  conclusions  are  derived  from  the  results  presented  in  Chapter  5.  The 
results  exhibit  new  ways,  not  predicted  by  analyses  using  steady  shears,  in  which  short  waves 
are  focussed  to  unstable  amplitudes.  These  new  ways,  grouped  with  the  critical-layer  interac¬ 
tion,  are  referred  to  as  refractive  convergence. 


CHAPTER  2 


BACKGROUND 


Two  mechanisms  proposed  for  the  dissipation  of  short  waves  involve  refraction  by 
steady  shears:  the  critical-layer  interaction  and  the  shear  instability  of  Phillips  (1966)  and  Hart¬ 
man  (197S).  These  mechanisms  are  examined  in  this  chapter. 

The  assumptions  of  inviscid,  small  amplitude  disturbances  in  a  steady,  horizontal 
current  yield  a  simple  model  of  a  critical  layer  that  is  described  by  the  Taylor-Goldstein  equa¬ 
tion.  Writing  u' -  (u’,0,w')  with  w'(x,z,t)  -  w(z)eakx~*,)  as  the  perturbation  velocity,  k  as 
the  horizontal  wavenumber,  a>  as  the  Doppler-shifted  frequency,  N  as  the  Vaisala  frequency,  U 
as  the  mean  flow  velocity,  and  c  —  tn/k  as  the  Doppler-shifted  horizontal  phase  velocity,  the 
Taylor-Goldstein  equation  becomes 


hu  + 


N 2 

( c-U )2 


+ 


Li¬ 
te- U) 


0  . 


(2.1) 


The  critical  layer  appears  as  a  singular  point  in  the  equation  and,  provided  U:  is  nonzero  at  the 
critical-layer  depth,  zc,  a  branch  point  in  the  general  solution.  A  series  solution  to  (2.1)  reveals 
that 


w-  A  (z-zc)'/i+,<*'“,/4,'A  +  B{z—zcyh~aRi~x/*)h ,  (2.2) 

as  z  —  Zf,  where  Ri^N2/U2  is  the  local  Richardson  number  and  A  and  B  are  complex  con¬ 
stants.  In  the  unstratifled  case,  N  -  0  and  the  addition  of  a  small  amount  of  viscosity  removes 
the  singularity.  With  stratification,  both  density  and  momentum  become  unbounded  at  the  crit¬ 
ical  layer,  so  both  viscous  and  thermal  diffusion  must  be  included  to  remove  the  singularity. 
The  result  is  then  a  sixth-order  differential  equation  free  of  singular  points  for  all  values  of  c, 
which  must  be  examined  in  the  limit  of  vanishing  viscosity  and  diffusion.  For  those  who  prefer 
second-order  differential  equations  to  sixth-order  ones,  two  popular  alternatives  exist,  both  of 
which  move  the  singularity  off  the  real  axis  into  the  complex  c-plane.  Rayleigh  damping  added 
to  the  momentum  and  density  equations  is  one  way  to  accomplish  this  (LeBlond  and  Mysak, 
1978),  but  the  proper  method  is  to  set  up  an  initial  value  problem  using  Laplace  transforms 
(Miles,  1961).  The  resulting  equation  need  not  be  solved,  only  examined  to  locate  the  new 
position  of  the  singular  points. 

Once  the  branch  is  chosen,  reflection  and  transmission  coefficients  can  be  computed 
for  waves  encountering  a  critical  layer.  The  overreflection  of  internal  waves  was  first  discovered 
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in  this  way  (Jones,  1968).  The  problem  with  this  type  of  analysis  is  that  although  the  vertical 
velocity  decreases  to  zero  at  the  critical  layer,  virtually  nothing  else  does:  the  horizontal  pertur¬ 
bation  velocity,  the  vertical  and  horizontal  displacements,  the  perturbation  pressure  and  den¬ 
sity,  the  action  density,  the  intrinsic  energy  density  (see  (5.11),  and  the  wave  induced  shear  all 
become  infinite,  indicating  the  inconsistency  of  the  steady  state  linear  theory.  For  this  reason, 
BB  resort  to  an  initial  value  problem. 

Despite  the  inconsistency,  much  effort  has  gone  into  solving  the  Taylor-Goldstein 
equation  with  various  mean  flow  profiles,  density  profiles,  and  boundary  conditions  for  applica¬ 
tion  to  internal  wave  propagation  and  to  the  stability  of  disturbances  in  a  stratified  fluid  when 
Ri  <  { 4  In  stability  theory,  solutions  with  unbounded  behavior  at  the  critical  layer  represent 
neutral  disturbances  and  lie  on  curves  in  the  k-Ri  plane  called  singular  neutral  curves.  Miles 
(1963)  proves  that  contiguous  to  these  curves,  unstable  solutions  exist.  Since  their  phase 
speeds  have  nonzero  imaginary  part,  the  branch  point  disappears  for  real  values  of  U  and  the 
inviscid  solution  is  bounded  through  the  depth  where  U ,  and  the  real  pan  of  c  match  (Gold¬ 
stein,  1931). 

When  applying  the  Taylor-Goldstein  equation  to  internal  wave  propagation,  the  prem¬ 
ise  is  that  the  steady-state  results  wilt  match  those  of  a  nonsingular  analysis  except  for  the 
region  near  the  critical  layer.  BB  establish  that  this  is  so  in  their  initial  value  problem:  "The 
important  point  is  that,  long  before  [the  linear  theory  becomes  invalid],  the  flow  above  and 
below  the  critical  layer  has  settled  down  to  a  steady  state  with  a  small  decaying  oscillation  super¬ 
imposed  and  that  the  steady  state  is  the  same  as  that  observed  by  ignoring  the  details  of  the 
critical  layer,  but  matching  around  it  ...”  The  steady  solution  predicts  a  reduction  in  the  flux  of 
pseudomomentum  tTST  across  the  critical  layer  by  a  factor  of  exp[2ir(/?t— 1/4)1'4]. 

The  dynamics  of  critical-layer  interactions  are  best  understood  within  the  WKB 
approximation  (Bretherton,  1966).  As  a  group  of  waves  propagates  toward  a  critical  layer, 
refraction  causes  the  intrinsic  frequency  and  group  velocity  to  decrease.  The  total  action  of  the 
group  remains  constant,  but  because  it  is  crowded  into  a  smaller  volume,  the  action  density 
grows.  Though  not  conserved,  the  intrinsic  energy  density  £  is  similarly  amplified  through 
convergence  induced  by  the  currents. 

A  second  type  of  dissipation  mechanism  for  short  waves,  originally  proposed  in  Phil¬ 
lips  (1966,  first  edition  only),  is  more  thoroughly  explored  by  Hartman  (1975).  Hartman  stu¬ 
dies  a  wave  packet  in  a  current  of  constant  shear  under  restrictions  that  limit  dispersion  and 
geometrical  spreading  due  to  the  variations  of  group  velocity  of  the  Fourier  components.  The 
intrinsic  frequency  at'  (related  to  the  Doppler-shifted  frequency  at  by  at  -  at’  +  Uk  for  k  >  0) 
and  the  vertical  wavenumber  m  satisfy 
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^r-  «  1  ,  —  «  1  ,  (2.3a,b) 

w,  mc 

where  6o>'  and  6m  are  the  changes  over  a  distance  of  one  inverse  wavenumber,  and  the  sub¬ 
script  c  denotes  the  center  value  for  the  packet.  A  fixed  observer  sees  roughly  the  same  intrin¬ 
sic  frequency  everywhere  but  sees  a  wide  range  of  Doppler-shifted  frequencies  because  U 
varies  over  the  vertical  dimension  of  the  packet.  In  the  steady  solutions  to  the  Taylor- 
Goldstein  equation,  the  Doppler-shifted  frequency  remains  constant  while  the  intrinsic  fre¬ 
quency  varies  with  depth. 

The  long-time  packet  behavior  differs  from  the  solutions  obtained  by  BB.  In  addition 
to  w',  u'  and  E  vanish  in  the  limit  of  large  time.  While  the  vertical  wavenumber  and  wave- 
induced  shear  increase,  by  the  time  they  become  sufficiently  large  to  allow  instability,  E  is 
greatly  reduced. 

The  cause  of  the  difference  is  that,  in  this  packet  setup,  the  ray  paths  are  not  converg¬ 
ing  to  a  single  critical-layer  position,  or  even  a  critical-layer  region  that  is  narrow  relative  to  the 
vertical  dimension  of  the  packet.  Instead,  lines  of  constant  phase  flatten  uniformly  in  depth  as 
a*'  decreases.  Ignoring  packet  edge  effects,  which  do  not  alter  the  conclusions  here,  this  absence 
of  convergence  reduces  the  equation  for  the  action  density  A  (see  (5.8))  to  dA/dt  -  0  because 
there  is  only  a  minor  gradient  in  A  or  in  the  group  velocity  across  the  packet;  accordingly,  as  «' 
decreases,  so  must  E. 

A  similar  asymptotic  behavior  governs  the  decaying  modes  associated  with  the  sudden 
start  from  rest  in  the  BB  investigation.  Why  these  two  limits  should  agree  is  explained  physi¬ 
cally  by  BB  and  mathematically  by  Brown  &  Stewartson  (1980).  An  initial  perturbation  can  be 
divided  into  a  discrete  spectrum,  with  finite  energy  in  a  mode  of  zero  bandwidth,  and  a  continu¬ 
ous  spectrum.  Given  enough  time,  the  disturbances  in  the  continuous  spectrum  always  decay 
in  a  constant-shear  flow  because  each  individual  component  selects  its  own  critical  layer.  Any 
initial  convergence  ends  as  the  rays  (in  a  WKB  sense)  become  more  horizontal. 

In  Hartman's  work,  the  decay  is  monotonic  because  of  the  initial  conditions.  These 
conditions  might  arise  as  the  mean  flow  accelerates  from  zero  velocity  through  a  sequence  of 
constant-shear  profiles  (Thorpe,  1978).  If,  on  the  other  hand,  a  packet  propagates  from  a 
motionless  region  into  a  shear  layer,  both  6w7o>'c  and  Ww’(  are  small.  Since  the  spread  in 
critical  layers  is  also  small,  ray  paths  initially  approach  each  other  possibly  producing  large 
values  of  E  before  decaying.  It  is  therefore  not  just  a  continuous  range,  but  a  continuous  and 
broad  range  of  Doppler-shifted  frequencies  that  differentiates  Hartman's  theory  from  the 
critical-layer  theory  of  Bretherton  (1966)  and  from  those  theories  based  on  the  Taylor- 


Goldstein  equation. 


Although  Hartman's  analysis  applies  only  to  currents  with  constant  shear,  an  important 
implication  is  not  so  restricted:  the  initial  orientation  of  the  rays  can  be  as  important  as  the 
structure  of  the  currents  for  short  wave  dynamics. 

An  objective  of  Chapter  5  is  to  use  ray  theory  to  contrast  the  dissipation  mechanisms 
reviewed  in  this  chapter  with  the  ways  in  which  inertial  currents  refract  short  waves  to  unstable 
amplitudes.  We  also  return  to  the  subject  of  packet  bandwidth  in  Chapter  5,  Section  C.  The 
accuracy  of  ray  theory  is  first  discussed  and  tested  in  the  following  two  chapters  under  the  con¬ 
ditions  observed  in  the  deep  ocean. 


CHAPTER  3 


% 


SCALE  ESTIMATES  AND  WKB  ACCURACY 


The  analyses  of  Hartman  (1975)  and  BB  apply  without  restriction  to  slowly  varying 
mean  flows,  but  since  the  goal  of  this  thesis  is  to  incorporate  inertial  shears,  which  are  neither 
steady  nor  constant  in  z,  progress  hinges  more  heavily  on  WKB  techniques. 

In  this  chapter  and  the  next,  an  attempt  is  made  to  obtain  some  indication  of  WKB 
accuracy  by  confining  attention  to  steady  mean  flows.  The  results  of  Hartman  (1975)  show 
negligible  amounts  of  partial  reflection  and  excellent  agreement  with  WKB  predictions.  The 
estimates  of  vertical  wavenumber  agree  exactly,  since  the  equation 


9m  ,9U 
9t  9z 


(3.1) 


holds  whether  the  mean  flow  varies  slowly  or  not.  Action  conservation  also  applies  generally, 
but  the  action  equation  becomes  approximate  for  the  case  studied  by  Hartman  when  A  is 
replaced  by  £W  (Andrews  &  McIntyre,  1978). 

A  further  step  toward  validating  WKB  theory  in  the  present  context  is  to  show  that 
partial  reflection  is  small  when  the  dynamics  are  described  by  the  Taylor-Goldstein  equation. 
Wave  properties  then  depend  strongly  on  z  for  typical  oceanic  conditions,  as  the  following  cal¬ 
culations  illustrate. 

If  the  shear  is  constant,  the  depth  of  the  turning  point  z,  is  related  to  the  depth  of  the 
critical  layer  zc  by 


j.  _  Oi'h  — 

rc  Ri  2n  , 


(3.2) 


where  A*  —  2 ir/k.  Equation  (3.2)  is  derived  using  o>  -  &/  +  Uk,  which  may  be  rewritten  as 


k(zc-z) 
S  “  Ri* 


(3.3) 


For  Ri  <:  10,  a  distance  of  less  than  half  a  horizontal  wavelength  separates  the  turning  point 
from  the  critical  layer  (in  contrast  with,  for  example,  LeBlond  and  Mysak,  1978,  figure  41.1). 

Another  measure  of  proximity,  which  is  more  pertinent  to  WKB  accuracy,  is  the  time 
required  for  a  wave  group  to  reach  the  critical  layer.  Within  the  WKB  approximation  this  time 
is  infinite,  but  by  integrating  along  a  ray  to  a  level  just  below  the  critical  layer  where  the 
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amplitudes  and  shears  are  large,  a  time  scale  shorter  than  a  few  Vaisala  periods  is  predicted  and 
confirmed  by  experiment  (Thorpe,  1981;  Koop,  1981). 

When  weighted  by  the  local  value  of  the  dimensionless  propagation  time  T. 
becomes  (for  constant  Uz ) 


-  /?/*log 


Ri* 


k(zc-z)  k2{zt-z)2 


Ri 


-  1 


(3.4) 


where  Ct j  is  the  vertical  component  of  the  group  velocity.  A  conservative  choice  of 
kiZf—z)  —  .1  (see  Thorpe,  1981,  figure  2)  and  Ri  —  4  yields  the  results  w'/yv-.O 5  (from  (3.3)) 
and  r.=7. 

The  effects  of  such  rapid  variation  on  WKB  accuracy  can  be  estimated  by  introducing 
the  variable  (CMver,  1974) 


£  —  jm(z')dz',  (3.5) 

where  m  is  the  coefficient  appearing  in  the  Taylor-Goldstein  equation,  (2.1).  and  r0  is  an  arbi¬ 
trary  depth.  Substituting  f  into  the  Taylor-Goldstein  equation  leaves 

w,i  +  +  -ir  4r-w?  “  0  .  (3.6) 

”  m 2  bz  5 


The  WKB  method  attempts  a  solution  for  w  in  terms  of  an  asymptotic  expansion.  To 
leading  order,  w  is  given  by 


M) 


(  m(zo ) 

V  m(z) 


e'f 


(3.7) 


provided 


9m 

bz 


«  1  . 


(3.8) 


When  the  shear  is  constant,  or  the  term  Us/(c—U)  in  m  is  negligible,  this  becomes 


1 

Ri* 


1 
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l_£Li 

.V2 


«  !  , 


(3.9) 
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which  is  satisfied  whenever  Ri  is  large  except  near  turning  points.  A  similar  constraint  on  the 
Richardson  number  is  derived  by  Hartman  (1975). 

The  criterion  (3.8)  emphasizes  that  the  validity  of  (3.7)  depends  on  the  fractional 
change  in  m,  not  m  itself,  and  that  the  fractional  change  must  be  small  over  a  distance  of  Mm, 
not  a  vertical  wavelength.  It  can  also  be  shown  that  the  expansion  is  doubly  asymptotic:  the 
error  approaches  zero  as  n  — *  0  for  fixed  z,  or  as  z  —  r0  for  fixed  fi  (Olver,  1974).  Stated 
another  way,  the  solution  (3.7)  is  accurate  through  a  nondimensionalized  (by  m)  distance  of 
the  order  of  l/fx.  Since  an  infinite  number  of  vertical  wavelengths  separate  z0  from  the  critical 
layer,  the  WKB  approximation  must  eventually  break  down  (Bretherton,  1966). 

Judging  from  Sanford’s  measurements,  the  condition  (3.9)  may  be  violated  frequently, 
even  away  from  turning  points.  Because  the  solution  for  w  given  by  (3.7)  is  the  first  term  in 
an  asymptotic  expansion,  the  result  could  be  improved  by  working  out  higher  order  terms.  If 
these  terms  are  small  relative  to  the  first,  it  would  be  tempting  to  conclude  that  the  first  term 
alone  is  adequate;  this  reasoning  may  cause  troubles,  however,  in  dealing  with  reflection.  Since 
the  energy  reflected  by  smooth  changes  in  the  medium  is  exponentially  small  as  n  —  0,  it 
never  appears  at  any  order  in  the  expansion  (Meyer,  1980).  Higher  order  terms  represent  small 
changes  in  phase  and  amplitude  but  never  generate,  for  example,  downgoing  waves  from  upgo- 
ing  ones. 

The  process  of  comparing  higher  order  terms  to  check  the  validity  of  the  WKB  method 
therefore  has  limitations.  Only  as  n  —  0  is  the  error  in  an  asymptotic  expansion  the  order  of 
the  first  term  excluded.  Partial  reflection  may  account  for  the  difference  between  a  solution 
accurate  asymptotically  and  a  solution  accurate  for  finite  n- 


CHAPTER  4 


REFLECTION 


A.  Theory 

The  circumstances  in  which  a  smooth  velocity  profile  induces  partial  reflection  are 
examined  in  this  chapter  by  calculating  exact  solutions  to  the  Tayior-Goldstein  equation  that 
can  be  interpreted  as  traveling  waves  as  z  — *  ±°°.  We  choose  as  a  convenient  mean  flow 
profile 


U  -  i/,jtanh/3z  , 

for  which  the  Tayior-Goldstein  equation  is 


(4.1) 


+ 


N2  2&2U(l-U2)  _  2 

( c-U )2  c-U 


w  —  0  . 


(4.2) 


Introducing  the  dimensionless  quantities 


Uo 


Uo 


H2  ± 

P2U$  '  “  "  P 


and  transforming  the  dependent  variable  from  z  to  U-  turns  (4.2)  into  a  Fuchsian  equation 
with  four  singular  points: 


wu.u.  ~ 


2U. 

1  -U.2 


wu.  + 


ZU . 


(c.-i/.)2d-C/.2)2  (c.-U.)d-U.2)  (l-U-2)2 


«•  0 


(4.3) 


The  singularities  are  located  at  {/.  —  —1,  1,  c.,  and  ».  Any  equation  with  four  regular  singular 
points  can  be  transformed  to  Heun’s  equation  by  a  two  step  procedure.  The  first  step  is  to 
bring  the  singularity  at  U.  -  -1  to  U-  ■  0  so  that  three  of  the  singularities  are  at  0,  1,  and  <». 
This  is  achieved  through  the  transformation 


U\  * 


(£M-1) 


Defining 


12 


13 


c  i 


(c.+l) 

2 


equation  (4.3)  becomes 


*'£/,*/,  + 


J 

16ci 


+ 


1 

y  a2(l-c,) 

1 

J 

—  2£/,  +  1 

£/.-! 

16(1 — C])  4 

(f/|-C,) 

16ci(l-C|) 

(4.4) 


The  solution  to  (4.4)  is  represented  by  the  Riemann  symbol 


0  1  Ci  oo 

w  <■  Z1!  i/i j  /7*2  i/4(1+/v)  2  £/ 1 

— i(i\  —ini  A(l—iv)  —1 


(4.5) 


where 


Ml 


M2  1 


2  1 4ci2 

1 

[  / 

2 

|  4(1— cj)2 

J 

4cf  (1— Cj)2 


—  a 2 


-  1 


(4.5a) 


(4.5b) 


(4.5c) 


The  exponents  of  (4.5)  govern  the  behavior  of  the  solution  at  each  of  the  singular 
points.  For  example,  if  A  and  B  are  complex  constants,  then 

w  ~  Alfr  + 


-  Aei<l'XnU'  +  Se"'1*’1”47' 


as  U\  —  0.  In  this  limit 


Milnt/i  -  mi  , 


so  fi\  is  identified  as  the  dimensionless  vertical  wavenumber  as  z  —  — Similarly,  n2 
corresponds  to  the  dimensionless  vertical  wavenumber  as  z  —  +<». 
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The  second  step  in  putting  (4.3)  in  the  form  of  Heun’s  equation  requires  a  transforma¬ 
tion  of  the  dependent  variable  w  by 

w,  -  Ur'Wrir'HUrcO-*""''*  .  (4.7) 


Inserting  this  into  (4.4)  leaves 

y  v+2^i+2/xj+3 


•*1 


u,u 


i"i 


Ux 


Ux-\ 


U\~cx 


5/3  U\-b 


t't  t/,(t/,-l)(C/,-c.) 


wi  -  0  , 


where 


-  5  ,  .  .  .  .  iv 

a  -  —  +  im  +  tn2  +  ~2  ' 

0  -  -Vi  +  ifi  i  +  />2  +  -y-  . 

y  «  1  +  2/>i  , 

8  ■»  1  +  iv 


and 


(4.8) 

(4.9) 

(4.10) 

(4.11) 

(4.12) 


b  «  l  +  +  ciOi^+Mf—  iipn+m)) 

-  Wl+M(l+/pi+Vi(l-/*’)(ci-l))  .  (4.13) 


(In  defining  the  parameters  in  (4.9)  to  (4.13),  The  notation  of  Snow  (1942)  is  followed  as 
closely  as  possible.)  The  Riemann  symbol  is  now 


1  C(  oo 

0  0  a  Ui 


(-2/**!  -2/>2  -iv  0 


(4.14) 


Equation  (4.8)  is  Heun’s  equation.  That  solution  which  is  analytic  at  U i-0  (i.e.  has 
zero  exponent)  is  known  as  Heun’s  function  and  is  given  by  the  series 


F(ci,b\a,0,y,SM\)  -  'Z,anU\  . 

0 


(4.15) 


A  three-term  recursion  formula  determines  the  a„'s: 

.  -b 

0o  "  1  .  ®i  -  — * 

yc\ 
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and 


0„+2(n+l-b')ci  - 


0«+t 


(/i+l)2(l+Ci)  +  (/i4*l)  l+(a+/3— 8)C|—fr)|  —  fl„(/f+a)(/i+0) 


(4.16) 


Drazin  (1958)  derives  an  equation  similar  to  (4.8)  with  c  -  0  to  study  shear  instability 
(see  equation  (23)  of  Drazin's  paper).  Drazin  notes  that  if  w,  is  constant,  a  solution  that 
satisfies  the  boundary  condition  of  wi  —  0  as  z  —  ±«>  is  obtained  if 

Sj§i/,  -  6  -  0  .  (4.17) 

In  Drazin’s  problem,  b  is  zero  and  a  can  never  vanish,  so,  by  equating  j§  to  zero,  the  singular 
neutral  curve 

J  -  a2(l-a2)  (4.18) 


is  found. 

Here  we  seek  solutions  to  (4.8)  that  represent  an  incident  plus  a  reflected  wave  as 
z - oq  (f/,-0) 

w(r)  -  e""*1*  +  Re'M|*  ,  (4.19) 


% 


and  a  transmitted  wave  as  z  — *  +<»  ({/,  —  1) 


w(z)  -  Te-^‘  . 


(4.20) 


The  amplitude  of  the  incident  wave  in  (4.19)  has  been  set  equal  to  unity,  and  the  complex  con¬ 
stants  R  and  T  are  related  to  the  fractional  reflected  and  transmitted  energy  densities  by 


Er  -  I*  I2 


Et  -  I  rl2 


(a2  +  M22) 
(a2  +  Mi2)  ' 


(4.21) 


(4.22) 
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The  analytical  continuation  formulae  necessary  to  determine  R  and  T  are  given  in 
Appendix  A;  however,  the  resulting  expressions  (a4)  and  (a5)  in  terms  of  Heun's  func¬ 
tions  are  not  very  useful.  A  simpler  and  more  interpretable  form  can  be  derived  if  the  number 
of  singularities  in  (4.3)  can  be  decreased  by  one,  because  Heun's  equation  then  reduces  to  a 
hypergeometric  equation.  One  way  to  accomplish  this,  valid  in  the  limit  of  J  —  °°  or  c  1  Uo, 
is  to  ignore  the  U-  term  in  the  Taylor-Goldstein  equation,  which  is  responsible  for  the  singu¬ 
larity  as  U'  —*  <*>.  Under  this  assumption,  the  equation  for  I/? I  derived  by  van  Duin  &  Kelder 
(1982)  and  converted  to  the  present  notation  is 

cosh2(-^0  +  sinh2ir(/x2~Mi) 

I A  I2  - - - .  (4.23) 

cosh2(-^)  +  sinh2ir(/42+/X|) 


B.  Discussion 

Two  limiting  forms  of  (4.23)  follow  easily.  The  first  is  the  limit  of  large  Richardson 
number.  If  this  is  achieved  by  letting  /V  —  »  or  £  —  0,  then  /*,,  /x2,  and  v  tend  to  infinity, 
leaving  the  approximate  solution 

|  A I  —  e  2**'  (4.24) 

(van  Duin  A  Kelder,  1982).  If  instead,  the  Richardson  number  is  increased  by  letting  Uo  —  0, 
(4.24)  still  holds  provided  a  short-wavelength  approximation,  2nn\  »  1,  is  made.  This  case 
illustrates  the  exponentially  small  magnitude  of  the  reflected  energy  in  the  WKB  limit,  as  men¬ 
tioned  in  Chapter  3. 

The  second  limit  of  interest  is  c  J  Uo-  By  again  imposing  2nm  »  1,  the  relationship 
(4.24)  is  obtained.  Partial  reflection  is  therefore  unimportant  for  waves  that  just  miss  satisfying 
critical-layer  conditions.  Mied  A  Dugan  (1975)  reach  the  same  conclusion  from  a  numerical 
integration  of  the  Taylor-Goldstein  equation. 

When  a  critical  layer  exists  within  the  shear  layer,  partial  reflection  is  again  unimpor¬ 
tant,  except  at  low  Richardson  number,  where  waves  may  be  overreflected.  BB  attribute  the 
reflected  energy  appearing  in  their  results  not  to  partial  reflection,  but  to  reflection  from  the 
discontinuity  in  their  velocity  prolile.  This  energy  is  derived  from  waves  that  pass  through  the 
critical  layer  twice,  once  before  and  once  after  reflection.  Each  time  the  energy  flux  is  dimin¬ 
ished  by  an  exponential  factor  (the  BB  transmission  coefficient),  yet  the  discontinuity  still 
accounts  for  the  major  portion  of  reflected  energy. 
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It  is  of  interest  to  see  whether  a  symmetric  profile  produces  the  same  behavior.  One 
such  profile  that  approaches  constant  values  as  r  — ■  is  U  -  Uosech2(flz).  The  resulting 
Taylor-Goldstein  equation  can  be  transformed  into  a  Fuchsian  one  by  the  same  change  of  vari¬ 
able,  V  -  tanh(j3z),  that  is  used  for  the  tanh  profile.  Here,  though,  the  number  of  singularities 
increases  to  five  (except  when  c  -  £/0  which  leads  to  one  irregular  and  three  regular  singular 
points.)  Rather  than  proceeding  by  the  method  outlined  in  Appendix  B,  it  is  simpler  to  model 
the  entire  coefficient  of  the  Taylor-Goldstein  equation 


m 


2 


N 2 

C c-U )2  (c-U) 


—  k2 


(4.25a) 


by 


m2  -  A  +  Bsech2(@z) 


(4.25b) 


with  A,B  >  0  (Epstein,  1930).  Neglecting  the  term  U-/(c-U )  in  (4.25a)  we  can  solve  for  U 
to  obtain 


U-c 


_ N _ 

(k2  +  A  +  Bsech2Bz)'h 


(4.26) 


The  form  of  the  reflection  coefficient  depends  on  the  nature  and  position  of  the  zeros 
and  singularities  of  m2  (Meyer,  1980).  For  the  model  (4.25)  to  be  appropriate,  the  zeros  must 
resemble  classical  turning  points  and  the  singularities  must  resemble  critical  layers.  The 
coefficient  m  vanishes  at 


z 


(—/it) 

2 


±  log( 


2B+A 


) 


(4.27) 


while  the  singularities  of  m  lie  at  z  —  ±itr/fi.  A  Taylor  expansion  verifies  that  the  proper  lim¬ 
its  result:  m2  has  simple  zeros,  and  its  only  singularities  are  second  order  poles.  (Note  that  the 
antisymmetric  model  of  Epstein  (1930),  m2  -  A  +  2?tanh(/3z)  gives  a  first  order  pole  at  the 
singularity  and  is  therefore  inappropriate  for  a  stratified  shear  flow;  it  is  pertinent,  however,  to 
unstratified  problems,  including  Rossby  wave  -  mean  flow  interactions.) 

Epstein  obtains 

1/?|2  — - s—(,r</- - —  (4.28) 

sin2(7rrf)  +  sinhJ(ir-4r) 
p 

for  the  reflection  coefficient,  where 


d  -  yi  -  '/4(1  +  BIB2)'* 


(4.29) 


18 


and 

a  -  AVl 3  .  (4.30) 

For  the  three  cases 

(0  0-0, 

(ii)  0  —  0,  irA'Vfi  »  1,  and 

(iii)  0  —  oo 

IRJ  takes  the  approximate  form 

|0| -sin(ir0)e-2'^  .  (4.31) 

The  first  two  limits  correspond  to  Ri  —  °°  and  the  last  limit  to  c  -  £/o-  As  for  the  case 
U  -  t/otanh(0r),  the  reflected  energy  is  exponentially  small  for  waves  that  have  phase  speeds 
barely  exceeding  the  mean  flow  maximum. 

A  final  remark  concerns  the  validity  of  the  WKB  approximation.  Although  partial 
reflection  does  not  impair  the  accuracy  of  WKB  solutions  in  the  limit  of  J  —  °°  or  cl  Uq,  WKB 
approximations  fail  when  applied  to  a  curved  velocity  profile  with  low  Richardson  number. 
Then  the  term  Us/(c-U)  dominates  the  others  in  the  coefficient  of  the  Taylor-Goldstein  equa¬ 
tion.  An  example  depicting  the  importance  of  this  term  for  a  tanh  velocity  profile  is  plotted  in 
figure  4.1.  Within  the  WKB  approximation,  m2  increases  monotonically  with  depth  because  the 
shear  is  always  positive,  but  what  the  approximation  predicts  to  be  refraction  toward  low  fre¬ 
quency  turns  out  to  be  refraction  toward  a  turning  point,  since  the  exact  value  of  m 1  becomes 
negative  where  the  curvature  term  gains  influence.  It  appears  that  this  concern  and  not  partial 
reflection  is  the  major  drawback  in  applying  WKB  techniques  to  flows  described  by  the  Taylor- 
Goldstein  model. 

However,  the  application  of  the  Taylor-Goldstein  model  to  short  waves  in  inertial 
currents  has  severe  limitations,  since  the  approximation  of  steady  shears  is  valid  only  for  times 
that  are  short  relative  to  an  inertial  period.  (Note  that  in  the  results  of  figure  1.1,  the  time 
dependence  of  the  currents  becomes  important  before  one-tenth  of  an  inertial  period  has 
elapsed.)  This  restricts  the  model’s  utility  to  waves  of  high  intrinsic  frequency.  Another  factor 
enters  though:  the  vertical  gradient  in  the  Doppler-shifted  frequency  u.  If  this  gradient  is  large, 
as  in  the  problem  of  Hartman  (1975),  the  Taylor-Goldstein  model  is  inappropriate  even  for 
short  times. 

Nevertheless,  the  Taylor-Goldstein  equation  is  useful  in  testing  WKB  accuracy.  It  pro¬ 
vides  the  Richardson  number  criterion  given  in  (3.9)  and  an  indication  of  the  importance,  or 
unimportance,  of  partial  reflection.  Such  results  are  difficult  to  obtain  from  analyses  that  include 
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Figure  4.1.  The  coefficient  of  the  Taylor-Goldstein  equation.  The  asymptotic  predicitions 
(dashed  lines),  obtained  by  neglecting  the  curvature  term  Ua/(c-U)  in  the  Taylor-Goldstein 
equation,  are  compared  with  the  exact  values  (bold  lines).  The  calculations  are  made  for 
7-2.0,  c/Ug-3.0,  («'/A0_«,-.3,  and  Ar//3— .  1.  Asymptotic  theory  does  not  anticipate  the  turn¬ 
ing  point  just  above  fiz  —0. 
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time-varying  shears,  although  some  progress  along  these  lines  can  be  made  by  extending 
Hartman’s  theory  to  a  time-dependent  current  of  the  form  U.  -  U.(t  only)  (Thorpe,  1978). 


CHAPTER  5 


SHORT-WAVE  PROPAGATION  IN  AN  INERTIAL  CURRENT 

Large-amplitude  inertial  oscillations  with  vertical  wavelengths  of  a  few  hundred  meters 
or  less  are  found  wherever  measurements  are  made  in  the  ocean.  They  constitute  the  mqjor 
shear  throughout  the  water  column,  strongly  refracting  short  waves  and  creating  the  possibility 
for  short-wave - critical-layer  interactions. 

The  effects  of  inertial  currents  on  the  propagation  of  short  waves  and  on  the  evolution 
of  critical-layer  encounters  are  analyzed  here  using  ray  theory.  Strictly  speaking,  the  introduc¬ 
tion  of  time-varying  shears  eliminates  singularities  associated  with  critical  layers;  at  issue,  how¬ 
ever,  is  whether  behavior  resembling  the  BB  predictions  can  still  be  achieved  and,  if  so,  which 
conditions  favor  its  occurrence. 

Rotation,  which  has  previously  been  ignored,  is  introduced  in  this  chapter.  Jones 
(1967)  analyzes  the  critical-layer  interaction  by  modifying  the  Taylor-Goldstein  equation  to 
include  rotation.  He  finds  that  rotation  introduces  two  additional  singularities  at  the  depths  at 
which  «'  —  ±/,  where  <*>'  and  /  are  the  intrinsic  and  inertial  frequencies  respectively.  He  also 
finds  an  exponential  decrease  in  the  vertical  flux  of  horizontal  pseudomomentum  (Andrews  & 
McIntyre,  1978)  across  the  singularity  at  at'  -  /.  However,  the  result  most  relevant  to  the 
present  purposes  is  that,  in  a  routing  fluid,  critical-layer  convergence  and  the  shear  insubility 
predicted  by  Hartman  (1975)  occur  as  «'  — ■  /.  (In  Hartman's  theory,  this  is  a  WKB  result, 
since  at'  is  meaningfully  defined  in  his  analysis  only  in  the  WKB  limit.) 

We  take  as  a  simple  model  of  an  inertial  current  the  standing  wave 

U  «'((/,  V,0)  -  (/osin/3r  (sin ft ,cos /r ,0)  (5.1) 

with  wavelength  -  2w/0  and  period  T,  —  2i r/f.  We  consider  short  waves  of  wavenumber 
K  -  (k,0,m)  and  assume  that  the  short-wave  amplitudes  are  independent  of  horizontal  posi¬ 
tion.  Then  the  y— component  of  the  inertial  velocity  V  does  not  enter  into  the  ray  equations 
and  therefore  need  not  be  considered  in  the  following  calculations. 

The  short  waves  have  intrinsic  frequency  at',  Doppler-shifted  frequency  at,  and  since  k 
is  assumed  positive, 

a*  -  a>'  +  Uk  .  (5.2) 

The  group  velocity  of  the  short  waves  is  C,  -  (Cfi,0,Cj3),  with 
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^  ,/  (N2-/2)  k2m 

CM3  -  TT~  ’ 


<u 


(5.3a) 

(5.3b) 


and  as  in  previous  chapters,  c'  and  c  are  the  intrinsic  and  Doppler-shifted  horizontal  phase 
speeds  respectively. 


The  minimum  Richardson  number  of  the  inertial  current,  J,  and  N/f  are  set  to 


J  - 


N 2 

put  “ 


(5.4a,b) 


The  magnitude  of  J  represents  a  compromise  between  the  measured  values  of  Sanford  (197$, 
figure  12),  which  fall  within  1  and  4  ninety  percent  of  the  time,  and  concerns  for  WKB  accu¬ 
racy  (see  (3.9)).  The  value  for  N/f  is  appropriate  to  a  depth  of  roughly  1000  meters  at  midla¬ 
titudes. 


A.  Ray  Theory 

Writing  the  dispersion  relation  as 

w(r,r)  -  n(K,z,r) 


PN2+ m2/2 
P+m2 


+  kU 


to  be  integrated  along  rays  defined 


are  (Garrett,  1968) 


dw  .3 U 

dt  dt 


dA 

dt 


-Ahc« 


+  n, 


(5.5) 


(5.6a,b) 


(5.7) 


(5.8) 


where  A  is  the  wave  action  density.  We  have  adopted  the  notation  of  Hayes  (1970),  using  0 
and  denoting  partial  derivatives  by  subscripts  when  the  independent  variables  include  K  as  well 
as  z  and  t.  The  variable  x  enters  only  in  the  equation  for  the  x -component  of  the  ray  path. 
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In  all  other  equations,  d/dt  can  be  reduced  to  d/dr  4-  C|3d/dz. 

Rather  than  integrating  along  adjacent  rays  to  determine  dm/d2  in  (S.8),  Hayes  sug¬ 
gests  advecting  the  value  of  dm/dz  along  the  single  ray  defined  by  (5.6).  This  introduces 
another  equation 

<s” 

(see  Hayes  (1970),  equation  (19)).  From  the  dispersion  relation,  we  find  that 


„-0, 

(5.10a) 

,0  2U 
dz2  ’ 

(5.10b) 

3/2)(l  -  m2Ob~'  +  iV2a~")  , 

(5.10c) 

where  in  (5. 10c) 


a  —  k2N2  +  m2/2 
b-k2  +  m2  . 

From  A ,  we  compute  the  intrinsic  energy  density  £,  which  is  given  by 


£-  y  (u'2  +  v'2  +  w'2  +  9'2/N2: 


•A  «' 


(5.10d) 

(5.10el 

(5.11) 

(5.12) 


and  which  satisfies  the  ray  equation 


+  (5.i3) 

dt  dz  * ut  dt 

It  marks  the  energy  density  measured  by  an  observer  moving  with  the  local  velocity  of  the 

mean  current. 

We  also  compute  the  pseudoenergy  density  e,  related  to  A  and  £  in  the  WKB  limit  by 

e  ■  «  A  -  «  — '  (5.14) 

CD 

The  ray  equation  for  e 

,515) 


i  mm**#’ 


shows  that  e  is  conserved  whenever  w  is  constant  along  the  ray,  or  equivalently,  whenever  the 
medium  is  time  independent. 


In  the  results  presented  in  this  chapter,  E  is  emphasized  over  the  quantities  A  and  e 
because  it  provides  the  best  indication  of  whether  or  not  the  waves  reach  large  amplitudes.  In 
this  respect,  A  has  the  disadvantage  of  reflecting  changes  in  &>'  in  addition  to  amplitude.  In  the 
case  studied  by  Hartman  (1975),  for  example,  A  remains  constant  as  the  wave  amplitudes  grow 
or  decay.  A  similar  concern  holds  for  e,  which  may  even  be  zero  or  negative. 

To  calculate  the  exact  Richardson  number  as  a  function  of  time  and  depth,  the  ampli¬ 
tude  and  phase  of  the  short  waves  must  be  specified.  If  u'  and  v'  are  the  horizontal  perturba¬ 
tion  velocities  and  p'  the  perturbation  density,  the  Richardson  number  can  be  expressed  as 


Ri 


PO  02 


-£-(£/+«'))  +  -f- 

dz  dz 


(Y+v')\ 


(5.16) 


Letting  u'  —  ugcos*,  (5.16)  becomes 


Ri  - 


(cos/Sz)2 


1  +  ^-ycos* 


1  +  2y(cos*sm//+-4sin*cos/f)  +  y2(sin2*+-^r-cos34) 

U  (I)  1 


(5.17) 


where 


uo  c'  jn  1 
c'  Uo  P  cos pz 


(5.18) 


The  quantity,  ujc\  is  set  to  .2  at  the  initial  coordinates  and  computed  at  later  times  from  the 
WKB  relationship 


E- 


(5.19) 


with  E  obtained  from  the  solution  for  wave  action,  (5.8).  The  minimum  and  maximum 
Richardson  numbers  are  determined  by  running  through  50  values  of  *  at  each  time-step. 

Equations  (5.6)- (5.9)  are  integrated  using  the  variable-order  method  of  Butirsch  and 
Stoer  (1966).  As  a  check,  the  formula  for  the  intrinsic  frequency  is  included: 


(5.20) 
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The  integration  ends  if  £  exceeds  50  times  its  initial  magnitude  or  if  Ri  <  1/4.  We 
shall  refer  to  this  process,  by  which  unstably  high  amplitudes  are  attained,  as  a  focussing  or 
convergence.  Once  started,  focussing  typically  proceeds  at  a  rapid  rate,  so  the  final  values  of 
at,  x,  and  t  are  seldom  sensitive  to  the  cutoff  value  for  £.  In  addition,  the  focussing  of  rays  is 
the  sole  cause  of  instability  in  all  cases  presented  in  this  chapter:  examples  of  shear 
intensification  in  which  the  rays  remain  parallel,  as  envisioned  by  Hartman  (1975)  and  Phillips 
(1966)  for  steady,  linear,  mean  flows,  were  never  found. 

B.  Results 

In  figure  5.1,  a  ray  calculation  begins  at  z  -  0,  t  -  0  with  a  Doppler-shifted,  horizontal 
phase  speed  of  .6i/<>  The  ray  path  would  head  asymptotically  toward  a  critical  layer  near 
z  -  ,1X,  if  the  mean  flow  was  steady  and  maintained  at  U  —  6'osin(/3z). 

Instead,  as  the  short  waves  propagate  upward,  the  inertial  wave  accelerates  causing 
c/U0  to  double  from  .6  to  about  1.2  (figure  5.1b).  At  the  same  time,  refraction  reduces  the 
intrinsic  frequency  to  less  than  2/  (figure  5.1c).  Despite  the  increase  in  vertical  wavenumber 
(figure  5. Id),  there  is  no  tendency  for  shear  instability  (figure  5.1j)  because  the  intrinsic  energy 
density  £  is  greatly  reduced  (figure  S.lh)  and  because  the  mean  shear  is  low  near  the  level  of 
one-quarter  inertial  wavelength. 

Since  internal  waves  of  low  frequency  and  high  wavenumber  propagate  with  small 
group  velocity,  vertical  progress  is  slowed.  Caught  in  the  region  of  maximum  currents,  they 
are  swept  backwards  in  the  decelerating  phase  of  the  inertial  oscillation  (figure  5.1a).  A  full 
inertial  period  elapses  without  any  sign  of  instability  or  wave  growth.  It  is  not  until  the  waves 
begin  to  emerge  from  the  region  of  maximum  currents  that  £  intensifies  (reaching  50  times  its 
initial  value  before  the  calculation  is  stopped);  however,  at  this  time,  the  intrinsic  frequency 
and  group  velocity  are  increasing  (figure  5.1c),  which  means  that  refraction  is  in  a  direction 
opposite  to  that  required  for  a  critical-layer  interaction.  The  strong,  time-dependent  shears 
have  therefore  focussed  an  accelerating  group  of  short  waves  to  unstable  amplitudes  at  a  depth 
where  critical-layer  conditions  are  not  satisfied. 

The  ability  of  inertial  currents  to  focus  short-wave  groups  that  are  accelerating  will 
appear  in  many  of  the  results  presented  in  this  chapter.  It  is  a  type  of  convergence  seldom  con¬ 
sidered  by  steady-shear  analyses.  The  reason  is  that  these  steady  analyses  generally  assume  or 
lead  to  solutions  in  which  the  Doppler-shifted  frequency  <u  is  independent  of  depth.  Examples 
include  the  study  of  waves  generated  by  flow  over  topography  or  by  shear  instability.  In  the 
steady-state  case,  with  o>  constant,  the  ray  equation  for  A  takes  the  simple  form, 
ACgi  —  constant ,  so  A  intensifies  only  as  the  wave  group  decelerates. 
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Figure  S.l.  Summary  of  a  ray-tracing  integration  for  short  waves  propagating  through  an  iner¬ 
tial  current  given  by  (5.1).  All  plots  have  the  same  vertical  axis.  At  the  initial  time,  r— 0, 
If— 6,  and  c/Un-.6  Time  is  marked  parametrically  in  (a)  by  dots  dong  the  ray  path  every 
quarter-inertial  period.  Minimum  and  maximum  Richardson  numbers  (computed  from  (5.17)) 
are  shown  in  (J)  by  dashed  and  bold  lines  respectively.  In  addition  to  the  graphs  referenced  in 
the  text,  the  following  are  included: 

(e)  m~2bm/bz,  provided  as  ah  indication  of  WKB  accuracy, 

(g)  the  vertical  fluxes  of  A  and  e,  and 

(/)  the  terms  contributing  to  changes  in  £  (see  (5.13)). 

Note  that  convergence  takes  nearly  Vh  inertial  periods  to  occur,  and  that  at  this  time,  u>'  and 
Cgj  are  increasing. 
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Figure  5.2.  Same  as  for  figure  5.1  except  that  a  steady  mean  current  has  been  added  to  the 
inertial  oscilUation.  The  total  background  velocity  in  the  x-direction  is:  U—  £/<£in/3zsin/r  +  rr , 
where  r2— .001  Af2. 
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In  the  example  pictured  in  figure  5.2,  ail  the  initial  conditions  are  the  same  as  in  figure 
5.1,  except  that  a  linear,  steady  mean  flow  of  the  form  Ut  -  rz  has  been  added  to  the  inertial- 
wave  background.  This  is  intended  to  model  the  shear  of  a  geostrophic  current  or  mesoscaie 
eddy.  The  parameter  r  is  chosen  such  that  the  Richardson  number  of  the  steady  flow  is  1000, 
but  the  results  are  similar  for  a  Richardson  number  of  100. 

With  this  modification,  the  waves  are  refracted  to  small  vertical  wavelengths,  but  there 
is  no  tendency  towards  shear  instability  such  as  that  which  develops  within  a  few  Vaisala 
periods  in  Hartman’s  results.  The  Richardson  number  never  gets  near  one-quarter,  in  part 
because  of  the  curvature  of  the  velocity  profile  and  in  part  because  BU/Bz  vanishes  twice  every 
inertial  period.  The  short  waves  are  carried  back  and  forth  by  the  currents  for  nearly  four  iner¬ 
tial  periods. 

Unlike  the  previous  two  cases,  short  waves  originating  near  the  depth  r  *■  . 5\,  find 
conditions  resembling  those  of  critical-layer  interactions  (figures  1.1,  5.3).  These  waves 
encounter  a  decelerating  mean  flow  with  positive  shear  (.BU/Bz  >  0),  similar  to  the  situation  in 
which  Thorpe  (1981)  generates  critical  layers. 

Despite  an  initial  short-wave  phase  speed  of  lU^  the  convergence  of  ray  paths  shown 
in  figure  5.3  develops  well  before  the  inertial  current  attains  its  maximum  velocity,  and  conver¬ 
gence  continues  at  depths  where  the  shear  becomes  negative  and  thus  where  critical-layer  con¬ 
ditions  cannot  be  met;  moreover,  for  the  calculations  summarized  in  figure  5.3,  the  value  of  c 
is  greater  than  2  U0  at  the  time  of  focussing  and  is  therefore  more  than  twice  the  local  value  of 
U. 

That  critical  layers  need  not  account  for  the  majority  of  convergences  predicted  by  ray 
theory  is  emphasized  by  figure  5.4,  which  summarizes  the  results  of  576  sets  of  ray-calculations. 
(The  symmetry  U(z,t)  -  U(z  +  A,/2 ,t  +  TJ2)  for  U  given  by  (5.1)  reduces  the  number  of 
independent  cases  by  a  factor  of  two.)  The  initial  conditions  cover  twelve  times  evenly  spaced 
over  an  inertial  period,  twelve  depths  evenly  spaced  over  an  inertial  wavelength,  and  four 
intrinsic  phase  speeds:  c'/Uq  -  .2,  .6,  1.0,  and  1.4.  The  initial  intrinsic  frequency  is  held  fixed 
at  at'  -  .2.  Each  integration  continues  for  two  inertial  periods  or  two  inertial  wavelengths 
unless  convergence  causes  the  intrinsic  energy  density  E  to  intensify  (see  the  final  paragraph  of 
section  5A).  If  this  occurs,  the  integration  stops  and  the  final  time  and  depth  relative  to  the 
phase  of  the  inertial  oscillation  are  marked  in  the  figure. 

The  shear  is  positive  at  the  times  and  depths  falling  in  the  hatched  sections  of  the  plot. 
Although  critical-layer  interactions  are  restricted  to  this  area,  to  classify  all  the  cases  in  the 
hatched  sections  as  critical-layer  interactions  is  to  include  behavior  that  differs  substantially 
from  the  results  of  BB  and  Thorpe  (1981):  sudden  increases  in  E  occur  at  intrinsic  frequencies 


Figure  S.4.  The  location  of  convergences.  Plus  marks  indicate  the  depth  relal.  /e  to  an  inertial 
wavelength  X,,  and  the  time  relative  relative  to  an  inertial  period  Th  at  which  the  ray  calcula¬ 
tions  predict  convergence  (F>5Q£(r-0)  or  Ri<  'A).  The  figure  displays  only  the  phase  of  the 
inertial  current  at  which  this  happens.  A  ray  which  reaches  f-1.5T  and  r— 1.75X,,  for  exam¬ 
ple,  is  given  the  coordinates  (.5,.7).  Cases  in  which  the  short  waves  propagate  for  longer  than 
two  inertial  periods  or  two  inertial  wavelengths  without  becoming  focussed  are  not  plotted.  The 
initial  values  of  c'/Uo  are  .2,  .6,  1.0,  and  1.4.  Of  576  cases,  434  points  are  plotted  (75%  of  the 
total)  with  246  points  falling  in  the  hatched  region,  where  critical-layer  interactions  are  possible. 
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much  greater  than  /,  at  times  when  the  current  is  weak  everywhere,  and  at  depths  where  the 
mean  shear  vanishes. 

The  number  of  cases  that  can  be  attributed  to  a  critical-layer  interaction  depends  on  its 
definition.  Let  a  critical-layer  interaction  be  one  in  which  wave  amplitudes  intensity  as  «'  —  / 
and,  in  the  WKB  approximation,  as  Ct j  —  0.  Then  only  10%  of  the  cases  plotted  in  figure  5.4 
fit  this  definition  if  by  &»'  —  /,  we  mean  that  «'  <  5/  when  focussing  occurs.  This  increases  to 
25%  if  the  condition  changes  to  ot  <  10/. 

Higher  percentages  follow  if  a  linear  mean  flow  is  added,  as  in  figure  5.2,  or  if  the  ini¬ 
tial  intrinsic  frequency  is  reduced;  however,  as  long  as  c/  Uq  <,  1  at  the  initial  position,  the 
merging  of  rays  continues  *  '  be  split  between  shears  of  positive  and  negative  sign. 

At  higher  phasr  speeds  (figure  5.5),  the  number  of  convergences  is  reduced.  Again, 
convergence  is  not  limited  to  regions  of  positive  shear  nor  to  waves  of  low  frequency.  The  final 
value  of  <■>’  falls  below  its  initial  value  of  10/  in  less  than  half  of  the  cases  plotted  in  figure  5.5. 
We  also  note  that  since  a  high  intrinsic  frequency  corresponds  to  a  low  ratio  of  m/k,  focussing 
affects  more  than  just  the  shortest  vertical  wavelengths.  If  focussing  develops  when  ot'  -  10/, 
for  instance,  it  takes  a  horizontal  wavelength  of  only  500  meters  to  give  a  vertical  wavelength 
of  100  meters. 

The  value  of  dm/dz  is  one  of  several  initial  conditions  left  unchanged  so  far.  It  has 
been  set  to  zero,  indicating  that  there  is  no  initial  convergence  to  a  bundle  of  rays.  To  test 
whether  the  results  depend  as  critically  on  this  parameter  as  they  do  in  the  case  of  steady 
shears,  where  the  value  of  dm/dz  accounts  for  the  difference  between  the  Hartman  and  BB 
solutions  (Chapter  2),  the  following  example  is  considered. 

Suppose  high  frequency  waves  can  propagate  a  vertical  distance  of  one  or  two  inertial 
wavelengths  before  the  current  changes  markedly.  If  the  waves  are  generated,  for  example,  by 
shear  instability,  then  «  is  constant  and  dm/dz  can  be  derived  from  the  dispersion  relation. 
For  short  times. 


dm  ^  dm  dot' 
dz  dot'  dz 


N1-/1  dU 
W1-/1)111  (N1  — a/2)  1/2  dz 


1/2  S^m 


(5.21) 


This  is  used  as  an  initial  condition;  for  longer  times,  o>  is  allowed  to  vary  and  the  system 
evolves  according  to  (5.6M5.9). 
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Figure  5.6  contrasts  this  initial  state  (figure  5.6a)  with  one  in  which  &>'  is  independent 
of  depth  (figure  5.6b).  The  integrations  begin  at  t  -  *A7J  with  oi/N  -  .2  and  c/£/0  -  1.4  in  the 
former  case  and  oi'/N  —  .2  and  c'/(/0  -  1.4  in  the  latter  case.  The  rays  originating  from 
.4  <  r  <  .75  in  figure  5.6a  illustrate  that  divergence  in  one  region  may  strengthen  convergence 
in  another.  Whether  the  amplitudes  increase  in  a  bundle  of  rays  is  thus  as  much  a  question  of 
where  the  rays  end  up  relative  to  the  phase  of  the  inertial  oscillation  as  the  degree  of  conver¬ 
gence  with  which  they  start. 

Integrations  based  on  (5.21)  and  run  at  12  starting  times  covering  an  inertial  period 
show  patterns  similar  to  those  plotted  in  figure  5.5.  In  part  this  is  because  the  differences 
between  the  constant-^  and  constant-w'  initial  conditions  decrease  as  the  magnitude  of  the  iner¬ 
tial  current  decreases;  however,  even  when  the  current  is  strong,  the  difference  in  the  location 
of  instabilities  relative  to  the  phase  of  the  current  tends  to  lessen  as  the  calculations  proceed 
beyond  an  inertial  period. 

The  sensitivity  of  the  ray  solutions  to  changes  in  Ar/y3  is  shown  in  figure  5.7,  which 
contains  the  ray  trajectories  for  the  interval  between  t  -  0  and  the  time  indicated  in  the  plots. 
The  value  of  k/p  labels  each  curve.  The  trajectories  for  k/p  -  .40  in  figures  5.7a,b  correspond 
to  those  in  figure  5.1a  and  in  figure  1.1  (with  z  -  .3\,  at  t  -  0)  respectively.  The  former  is  an 
example  in  which  convergence  is  not  associated  with  critical  layers,  while  the  latter  models  the 
critical-layer  interaction  as  closely  as  any  of  the  results  shown  thus  far. 

The  final  spread  of  trajectories  in  figure  5.7a  is  less  than  ha»f  a  vertical  wavelength  of 
the  short  waves  (approximately  .2X,  at  t  -  1.471),  but  the  solutions  are  more  sensitive  to  varia¬ 
tions  in  k/p  than  is  revealed  by  the  trajectories.  For  k/p  ■■  .44,  .40  convergence  develops  at 
r  -  1.47),  while  for  k/p  —  .38,  .36  E  normalized  by  its  initial  value  remains  small  for  over  2  Vi 
inertial  periods.  This  dramatic  change  in  the  intrinsic  energy  density  for  slight  changes  in  the 
independent  parameter  is  a  characteristic  feature  of  the  ray  solutions. 

In  figure  5.7b,  convergence  occurs  in  ail  three  cases  at  approximately  the  same  depth 
and  time  (f  -  .1 7))  despite  a  spread  of  ±25%  in  horizontal  wavenumber.  This  suggests  that 
convergence  in  this  example  is  not  limited  to  narrow-banded  wave  packets.  The  effects  of 
bandwidth  are  considered  in  the  next  section. 

C.  Dispersion  and  geometrical  spreading 

Dispersion  results  from  the  dependence  of  Cf  on  the  magnitude  of  K,  while  geometri¬ 
cal  spreading  results  from  the  dependence  of  C;  on  the  direction  of  K.  These  processes  are 
analyzed  in  this  section  using  ray  theory  to  determine  their  effects  of  the  convergence  of  local¬ 
ized,  short-wave  packets. 


r/9z  takes  the  value  obtained  in  (3.13),  corresponding  to  <u  independent  of  depth,  while 
dm/dz- 0,  corresponding  to  <■>'  independent  of  depth. 
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The  importance  of  dispersion  and  geometrical  spreading  depends  on  the  wavenumber 
bandwidth  of  the  packet.  Bandwidth  has  been  ignored  thus  far,  under  the  premise  (verified 
below)  that  the  results  capture  the  essential  features  of  at  least  narrow-banded  wave  packets. 
With  bandwidth  neglected,  changes  in  the  action  density  A  have  been  attributed  entirely  to 
refraction,  as  though  ail  the  energy  resides  in  a  single  value  of  K  and  <*>  at  each  depth  and  time. 
A  finite  bandwidth  is  implied,  however,  since  the  action  propagates  at  the  group  velocity,  not 
the  phase  velocity. 

Ray  theory  can  be  used  to  study  the  processes  of  dispersion  and  geometrical  spreading. 
Suppose,  for  example,  that  a  source  localized  in  space  and  time  generates  a  broad-banded 
wavefield,  as  in  the  "seas’  of  surface  waves.  As  the  Fourier  components  spread  due  to  their 
difference  in  group  velocities,  as  in  the  "swell"  region,  a  gradient  in  K  develops.  Strictly  speak¬ 
ing,  all  wavenumbers  are  present  everywhere,  but  to  a  first  approximation,  one  wavenumber  is 
dominant  at  each  position  with  all  others  destructively  interfering.  This  gives  rise  to  a  nonzero 
VK  even  if  the  medium  is  uniform,  which  when  substituted  into  the  ray  equation  for  the 
action  density  A  describes  the  changes  in  packet  shape  resulting  from  dispersion  and  geometri¬ 
cal  spreading.  Thus  bandwidth  information  enters  through  the  value  of  VK. 

This  discussion  is  exactly  the  physical  justification  for  the  method  of  stationary  phase, 
and,  in  fact,  ray  theory  can  be  used  to  derive  the  stationary-phase  result 

A  -  r"  ,  (5.22) 

for  large  r,  where  n  is  the  number  of  spatial  dimensions.  To  demonstrate  this,  we  consider  the 
simple  case  of  one-dimensional  propagation  in  a  uniform  medium.  (The  generalization  to  n 
dimensions  is  given  in  Appendix  B.)  Since  the  wavenumber  k  is  then  conserved, 

k  -  k(x—Cgt)  ,  (5.23) 

where  C(  is  the  group  velocity.  Implicit  differentiation  with  respect  to  x  leaves 


where  O k/bx)0  is  the  value  of  bk/bx  at  t  —  0.  Substituting  (5.24)  into  the  one-dimensional 
action  equation 


(see  (5.8))  gives  in  the  long-time  limit 
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dA 

dt 


t  ' 


(5.25b) 


and  leads  to  the  asymptotic  behavior  given  by  (5.22)  provided  n«,  Bk/Bx  >  0  at  r  —  0. 


The  advantage  that  ray  theory  has  over  the  method  of  stationary  phase  is  that  it  is 
easily  extended  to  nonuniform  media.  Complications  arise  because  in  general  Bk/Bx  is  not 
known  analytically  and  because  in  three  dimensions,  (5.24),  or  (5.10),  must  be  replaced  by  the 
system  of  equations 


<T7  K 
dt 


-vk  vk  -  vk  nlK  -  nKx  vk  -  n„  , 


(5.26) 


(Hayes,  1970).  (Only  six  of  the  equations  in  (5.26)  are  independent  since  the  tensor  VK  is 
symmetric.)  The  components  of  which  are  listed  at  the  end  of  Appendix  B,  include 
fl«,  fi and  so  forth,  and  similarly  for  fllK  and  ft„. 

We  use  (5.26),  the  action  equation  in  three  dimensions 

-  -A  V K  :  nKa  ,  (S.27) 

and  the  ray  equations  (5.5)-(5.7)  to  include  the  effects  of  dispersion  and  geometrical  spreading 
with  those  of  refraction. 

The  results  are  presented  in  figure  5.8,  which  shows  the  intrinsic  energy  density  E 
normalized  by  its  initial  value  as  a  function  of  time  for  several  ray  calculations.  The  number 
that  labels  each  curve,  which  for  later  reference  we  call  y,  has  the  same  magnitude  as  each  of 
the  terms 


JL_  dm  J_  |  dk  3k  dk  3/  d/ 1  ... 

m2  dr  ’  k2  I  dx  ’  By  ’  3r  ’  By  '  Bz  j 

and  has  positive  (negative)  sign  when  the  rays  are  initially  divergent  (convergent).  For  initial 
divergence  (convergence),  the  sign  of  each  component  of  VK  is  chosen  to  have  the  same 

(opposite)  sign  as  the  corresponding  component  of  ftKK,  so  that  all  nine  terms  contributing  to 

the  dyadic  product  in  (5.27)  have  the  same  sign.  This  maximizes  the  effects  of  dispersion  and 
geometrical  spreading. 

In  figure  5.8a,  the  initial  conditions  are  those  used  in  figure  5. 1  except  for  the  value  of 
y.  (The  bold  curve  corresponds  to  the  calculation  in  figure  5.1.)  For  y  >  0,  dispersion  and 
geometrical  spreading  are  important  for  times  less  than  about  half  an  inertial  period.  Refractive 
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Figure  5.8.  The  effects  of  dispersion  and  geometrical  spreading.  Each  curve  is  labeled  by  the 
value  of  y.  The  magnitude  of  y  is  the  magnitude  of  each  of  the  terms  in  (5.28),  while  the  sign 
of  y  is  positive  for  initially  divergent  rays  and  negative  for  initially  convergent  rays. 

The  ray  path  and  the  values  of  u>\  c,  CfJ,  and  m  for  (a)  are  shown  in  figure  5.1.  For 
(6),  the  ray  path  is  the  one  originating  at  z  -  .3X,  in  figure  1.1. 


convergence  becomes  the  dominant  process  after  1.27],  and  the  results  illustrate  the  explosive 
growth  in  E  that  this  convergence  produces.  As  in  nearly  ail  the  . ray  solutions  summarized  in 
this  chapter,  there  is  no  question  as  to  when  focussing  occurs. 

In  the  case  where  y  —  -.5  in  figure  5.8a,  the  rays  are  so  sharply  focussed  at  the  outset 
that  unstable  amplitudes  are  shortly  reached.  In  the  Taylor-Goldstein  model,  this  degree  of 
convergence  is  produced  by  a  current  with  low  Richardson  number  (see  (3.9)).  In  the  present 
context,  the  focussing  arises  because  the  initial  conditions  place  Fourier  components  with  slow 
group  velocities  ahead  of  those  with  fast  group  velocities.  The  ray  equations  then  predict  con¬ 
vergence  without  background  currents. 

The  second  example  (figure  5.8b)  shows  the  effects  of  dispersion  and  geometrical 
spreading  in  a  situation  resembling  the  critical-layer  interaction.  Here  the  initial  conditions 
(except  for  the  values  of  y)  are  the  same  as  for  the  ray  that  begins  at  z  -  .3\,  in  figure  1.1. 
The  wave  group  thus  encounters  a  decelerating  current  with  positive  shear,  as  in  the  experi¬ 
ments  of  Thorpe  (1981).  In  the  case  with  no  initial  convergence  (the  bold  curve  in  figure 
5.8b),  unstable  amplitudes  are  reached  after  only  one-tenth  of  an  inertial  period. 

For  y  -  ±.l  this  result  does  not  change  appreciably,  so  the  case  of  small  bandwidth  is 
similar  to  the  case  of  zero  bandwidth.  This  agrees  with  the  conclusions  of  BB,  who  claim 
(page  535)  that  "...a  wave  packet  clustered  around  a  single  wave-number  Ar0  does  not  behave 
qualitatively  differently  from  a  single  wave-number...*,  and  the  result  is  expected  from  the 
laboratory  study  of  Thorpe  (1981).  In  Thorpe’s  experiments,  a  rough  estimate  of  the 
bandwidth  in  horizontal  wavenumbers  Bk/k  is  l/nx  where  nx  is  the  number  of  horizontal 
wavelengths  to  the  corrugated  lower  boundary  of  Thorpe’s  tank.  Thorpe  uses  the  two  values 
«x  -  8  and  nx  -  16.  In  addition,  there  is  the  bandwidth  in  frequency  associated  with  the  time 
dependent  current  in  Thorpe's  study. 

For  the  larger  value  of  y  -  .5,  the  spreading  of  Fourier  components  dominates  the 
convergence  due  to  refraction  by  the  currents.  This  example  resembles  the  Han  man  result 
more  than  the  critical-layer  interaction  (for  the  times  shown  in  the  figure)  because  the  energy  is 
not  being  concentrated  at  one  depth.  The  intrinsic  energy  density  remains  small  relative  to  its 
initial  value  until  t  -  7],  when  convergence  occurs. 

As  mentioned  above,  the  basis  for  using  ray  theory  in  these  examples  is  the  assump¬ 
tion  that  the  wave  field  can  be  described  locally  by  one  value  of  K  and  <o  at  each  time;  how¬ 
ever,  measurements  suggest  that  the  deep-ocean  internal  wave  field  more  closely  resembles  the 
'seas'  of  surface  waves  more  than  the  'swell”  since  many  wavenumbers  and  frequencies  exist  at 
each  position.  We  must  therefore  deal  with  rays  that  start  at  the  same  location  with  different 
wavenumbers,  as  pictured  in  the  plots  of  figure  5.7,  and  with  rays  that  arrive  at  the  same 
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position  and  time  through  different  trajectories. 

The  quantitative  analysis  of  a  short-wave  spectrum  in  inertial  shears  will  not  be  con¬ 
sidered  here,  except  to  mention  that  in  treating  a  spectrum,  it  is  convenient  to  introduce  a  spec¬ 
tral  density  for  action,  which  we  denote  by  »(K,x,/)  to  distinguish  it  from  A  (x,r),  the  action 
density  of  a  single  wavetrain  in  physical  space.  The  equation  governing  conservation  of  n  is 
(Holloway  1981) 


n,  +  SikjnXj  -  SlXjnk.  -  0  (5.29) 

which  reduces  to  (5.8)  when 

n(K,x,/)  -  A  (x,r)8(K-K(x,r))  (5.30) 

(J.  Lemer,  unpublished  manuscript).  The  subscript  j  in  (5.29)  indicates  a  summation  over  the 
components  of  K  and  x. 

By  regarding  the  phases  of  the  short  waves  as  random,  the  action  and  energy  densities 
in  physical  space  may  be  written  as 

/  /»(K,x,/)rfK  (5.31) 

J*n(k,x,/)n'(k,x,r)</K  (5.32) 

respectively,  where  SI'  is  the  intrinsic  frequency.  It  is  the  value  of  the  integrals  in  (5.31)  and 
(5.32)  and  not  the  spectral  densities  by  themselves  that  should  determine  whether  unstable 
amplitudes  are  reached.  The  evaluation  of  these  integrals  would  seem  to  be  a  necessary  step 
toward  understanding  the  mechanisms  by  which  internal  waves  dissipate. 


CHAPTER  6 


CONCLUSION 

When  the  background  current  is  an  inertial  oscillation,  the  process  by  which  short 
waves  are  focussed  to  high  amplitudes  differs  in  significant  ways  from  the  steady-shear  case.  As 
the  ray  tracing  calculations  of  Chapter  S  emphasize,  a  critical-layer  interaction  offers  only  one 
example  of  convergence,  one  in  which  refraction  leads  to  a  decrease  in  the  intrinsic  frequency 
toward  /  and  the  group  velocity  toward  zero.  In  inertial  shears,  other  possibilities  are  equally 
likely.  The  intrinsic  frequency  may  be  high  or  increasing,  so  that  wave  groups  are  accelerating 
when  convergence  develops.  Inertial  shears  also  can  leave  the  rays  sharply  focussed  at  depths 
where  BU/dz  vanishes  and  at  times  when  U  is  identically  zero.  The  magnitude  and  sign  of  the 
shear,  therefore,  provide  little  indication  of  the  likeliness  for  short-wave  convergence  when 
c^U0  (figure  5.4). 

The  results  also  differ  from  those  of  steady-shear  analyses  in  that  the  condition 
c  <  Uo  does  not  guarantee  that  the  waves  will  shortly  attain  high  amplitudes  or  develop 
unstable  shears.  This  is  due  to  the  reversing  nature  of  the  inertial  currents,  both  in  space  and 
time,  which  can  delay  or  prevent  convergence  and  allow  the  short  waves  to  persist  for  several 
inertial  periods  without  becoming  unstable  (figure  5.2).  By  contrast,  instabilities  arise  within  a 
few  Vaisala  periods  in  the  studies  of  Fritts  (1978),  Hartman  (1975),  and  Thorpe  (1981).  Fritts 
(1982)  points  out  this  difference  between  steady  and  time  dependent  currents. 

The  reversing  currents  also  negate  Hartman’s  mechanism  of  dissipation  (shear  instabil¬ 
ity  without  an  accompanying  increase  in  £).  However,  short  waves  may  induce  instability 
without  reaching  large  amplitudes  if  they  propagate  into  an  inertial  current  that  is  by  itself  mar¬ 
ginally  stable. 

The  conditions  that  give  rise  to  critical-layer  interactions  match  those  studied  by 
Thorpe  (1981):  a  decelerating  current  coupled  with  a  mean  shear  of  the  sign  necessary  to 
decrease  the  intrinsic  frequency  of  the  short  waves  in  the  direction  of  propagation.  This  combi¬ 
nation  of  deceleration  and  refraction  leads  to  a  rapid  growth  of  the  short-wave  amplitudes. 
Here  again,  though,  differences  appear  in  relation  to  steady  shears,  since  wave  groups  with 
phase  speeds  exceeding  the  mean-flow  maximum  by  two  or  three  times,  and  the  local  mean- 
flow  velocity  by  much  more,  are  focussed  to  unstable  amplitudes  by  the  inertial  currents. 

These  results  indicate  the  difficulty  in  defining  a  critical-layer  interaction  quantitatively. 
The  criterion  c  -  U,  which  suffices  in  the  Taylor-Goldstein  model,  is  artificial  in  the  present 
context,  because  waves  of  higher  phase  speeds  are  similarly  amplified  by  the  currents.  There  is 
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no  obvious  limit  to  the  values  of  c  or  the  intrinsic  frequency  «'  at  the  time  of  focussing  that 
can  be  used  to  distinguish  the  critical-layer  interaction  from  other  convergences.  One  possibil¬ 
ity  is  to  associate  with  critical  layers  all  convergences  with  a  decreasing  m  along  the  ray  (i.e.  ail 
those  falling  in  the  hatched  sections  of  figures  S.4  and  5.5),  but  this  would  include  behavior 
much  different  from  the  results  of  BB  for  the  reasons  mentioned  above  and  in  Chapter  5.  To 
avoid  the  difficulties  in  definition,  and  to  designate  all  focussings  induced  by  currents  whether 
associated  with  critical  layers  or  not,  we  use  the  term  refractive  convergence. 

The  qualitative  features  of  refractive  convergence  are  shown  schematically  in  figure 

6.1.  The  arrows  represent  a  pair  of  adjacent  group  velocity  vectors,  whose  orientations  change 
through  refraction  by  an  inertial  current.  This  example  corresponds  to  the  situation  in  figure 

5.1,  although  the  range  of  angles  of  the  vectors  relative  to  the  vertical  and  relative  to  each 
other  is  exaggerated  for  illustration  purposes. 

The  vectors,  which  begin  parallel,  remain  so  in  the  second  diagram  of  figure  6.1 
(r  -  7’/8).  The  action  density  stays  approximately  constant  during  this  time  interval  (cf.  figure 
5.10,  so  the  dynamics  are  similar  to  those  predicted  by  Hartman  (1975),  not  by  the  Taylor- 
Goldstein  model.  Divergence  begins  thereafter  and  continues  when  the  current  vanishes  at 
r-7)  (again  unlike  the  Taylor-Goldstein  model).  Convergence  takes  over  after  t  -  %Th  inten¬ 
sifying  in  the  last  diagram  of  the  sequence  at  a  level  where  the  shear  is  negative. 

Once  underway,  refractive  convergence  typically  leads  to  sudden  increases  in  E - 

so  sudden,  in  fact,  that  a  dear  distinction  exists  between  those  waves  that  reach  high  ampli¬ 
tudes  and  those  that  do  not.  Because  there  are  few  intermediate  cases,  the  findings  are  not  sen¬ 
sitive  to  the  criterion  for  identifying  a  convergence  event  (page  25). 

The  calculations  of  Chapter  5  are  based  on  ray  theory.  Chapters  3  and  4  are  devoted  to 
assessing  the  errors  in  this  approach.  A  comparison  with  solutions  to  the  Taylor-Goldstein 
equation  finds  partial  reflection  to  be  unimportant  for  waves  with  a  Doppler-shifted  phase  speed 
barely  exceeding  the  mean  flow  maximum;  however,  WKB  techniques  perform  poorly  when  the 
curvature  term  UJic—U)  dominates  the  others  in  the  coeffirient  of  the  Taylor-Goldstein 
equation.  The  solutions  then  resemble  shear  modes  (solutions  to  the  Rayleigh  equation)  rather 
than  internal  wave  modes,  and  the  internal  wave  dispersion  relation  is  of  little  use,  even  quali¬ 
tatively.  The  example  shown  in  figure  (4.1)  illustrates  that  a  turning  point  can  appear  at  a 
depth  where  the  ray  equations  predict  a  small  and  decreasing  intrinsic  frequency. 

Chapter  4  finishes  with  a  discussion  of  the  inappropriateness  of  the  Taylor-Goldstein 
model  for  short  waves  in  a  inertial  current.  The  reasons  include  the  time  dependence  of  the 
shear  and  the  nonzero  vertical  gradient  in  the  Doppler-shifted  frequency  <v  that  a  time- 
dependent  shear  induces. 
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Figure  6.1.  The  process  of  refractive  convergence.  The  illustration  shows  schematically  the 
changing  orientation  of  a  pair  of  adjacent  group  velocity  vectors.  The  conditions  correspond  to 
those  depicted  in  figure  5.1.  The  lengths,  angles,  and  vertical  separation  of  the  vectors  are  not 
properly  scaled.  The  horizontal  separation  need  not  he  considered,  since  the  wave  amplitudes 
and  inertial  current  do  not  depend  on  horizontal  position.  The  profiles  of  U  are  indicated  by 
dashed  lines,  and  the  time  in  inertial  periods  is  given  below  each  diagram. 

The  vectors  begin  parallel,  diverge  for  times  between  7)/8  and  37)/4,  and  converge  for 
later  times.  Convergence  strengthens  in  the  final  diagram  at  a  depth  where  the  shear  is  nega¬ 
tive. 
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Near  convergence  zones,  the  equations  fail  because  of  the  small-amplitude  approxima¬ 
tion.  Numerical  solutions  to  the  nonlinear  equations  and  laboratory  experiments  show  over¬ 
turning  in  the  vicinity  of  the  critical  layer,  but  in  a  "noisy"  ocean  wave-wave  interactions  may 
limit  wave  growth  before  the  critical  layer  is  reached.  The  conditions  for  triad  instability,  for 
instance,  may  be  realized  before  those  for  gravitational  or  shear  instability.  Holloway  (1981, 
1982)  addresses  these  points  and  cautions  against  ignoring  the  nonlinear  interactions  of  a  ran¬ 
dom  wave  field  in  critical-layer  studies. 

There  are  also  linear  processes  to  be  examined  that  may  affect  the  focussing  of  short¬ 
wave  energy.  For  example,  it  should  be  useful  to  analyze  with  the  ray  calculations  a  broad 
banded,  short-wave  spectrum  propagating  through  random,  low-frequency  shears. 

The  results  presented  here  cover  a  very  limited  range  of  the  numerous  independent 
variables,  which  include  /,  0,  N/f ,  and  the  initial  values  of  a>\  m,  and  dm/d:.  Although 
critical-layer  interactions  similar  to  those  studied  by  Thorpe  (1981)  can  be  identified  in  inertial 
currents,  these  results  show  that  convergence  may  develop  in  ways  that  do  not  resemble 
critical-layer  interactions  at  all.  The  calculations  reveal  the  importance  of  the  more  general  pro¬ 
cess  of  refractive  convergence,  which  operates  under  a  broad  range  of  conditions.  Because  of  its 
connection  with  dissipation,  refractive  convergence  merits  attention  for  the  same  reasons  previ¬ 
ously  attached  to  critical  layers  by  BB,  Thorpe  (1981),  Hartman  (1975),  Fritts  (1978,  1982), 
Munk  (1980)  and  others:  to  improve  our  understanding  on  the  subjects  of  short-wave  dynam¬ 
ics,  the  internal-wave  energy  budget,  and  mixing  processes  in  the  deep-ocean. 


References 


Andrews,  D.G.  Sc  ME.  McIntyre  1978  On  wave-action  and  its  relatives.  J.  Fluid  Mech.  89, 
647-664. 

Booker,  J.R.  Sc  F.P.  Bretherton  1967  The  critical  layer  for  internal  gravity  waves  in  a  shear 
flow.  J.  Fluid  Mech.  27,  513-539. 

Bretherton,  F.P.  1966  The  propagation  of  groups  of  internal  gravity  waves  in  a  shear  flow. 
Quart.  J.  Roy.  Meteor.  Soc.  92,  466-480. 

Brown,  S.N.  Sc  K.  Stewartson  1980  On  the  algebraic  decay  of  disturbances  in  a  stratified  fluid. 
/  Fluid  Mech.  100,  811-816. 

Bulirsch,  R.  Sc  J.  Stoer  1966  Numerical  treatment  of  ordinary  differential  equations  by  extra¬ 
polation  methods.  Numerische  Math.  8,  1-13. 

Drazin,  P.G.  19S8  The  stability  of  a  shear  layer  in  an  unbounded  heterogeneous  inviscid  fluid. 
J.  Fluid  Mech.  A,  214-224. 

Epstein,  P.S.  1930  Reflection  of  waves  in  an  inhomogeneous  absorbing  media.  Proc.  Nat. 
Acad.  Set  16,  627-637. 

Fritts,  D.C.  1978  The  nonlinear  gravity  wave-critical  layer  interaction.  J.  Atmos.  Sci.  35,  397- 
413. 

Fritts,  D.G  1982  The  transient  critical-level  interaction  in  a  Boussinesq  fluid.  J.  Geophys.  Res. 
87,  7997-8016. 

Garrett,  GJ.R.  1968  On  the  interaction  between  internal  gravity  waves  and  a  shear  flow.  J. 
Fluid  Mech.  34,  711-720. 

Goldstein,  S.  1931  On  the  stability  of  superposed  streams  of  fluids  of  different  densities.  Proc. 
Roy.  Soc.  Load.  A132,  524-548. 

Hartman,  R.J.  1975  Wave  propagation  in  a  stratified  shear  flow.  J.  Fluid  Mech.  71,  89-104. 
Hayes,  W.D.  1970  Kinematic  wave  theory.  Proc.  Roy.  Soc.  Load.  A  320,  209-226. 

Holloway,  G.  1980  Oceanic  internal  waves  are  not  weak  waves.  J.  Phys.  Oceanog.  10,  906-914. 
Holloway,  G.  1981  Theoretical  approaches  to  interactions  among  internal  waves.  In  Nonlinear 


46 

Properties  of  Internal  Waves,  A.I.P.  Conference  Proceedings  No.  76,  American  Institute  of 
Physics,  New  York,  37-67. 

Jones,  W.L.  1967  Propagation  of  internal  gravity  waves  in  fluids  with  shear  and  rotation.  J. 

Fluid  Mech.  30,  439-448. 

Jones,  W.L.  1968  Reflexion  and  stability  of  waves  in  stably  stratified  fluids  with  shear:  a 
numerical  study.  J.  Fluid  Mech.  34,  609-624. 

Koop,  C.G.  1981  A  preliminary  investigation  of  the  interaction  of  internal  gravity  waves  with 
a  steady  shearing  motion.  J.  Fluid  Mech.  113,  347-386. 

LeBlond,  P.H.  &  L.A.  Mysak  1978  Waves  in  the  Ocean.  Elsevier  Scientific  Publishing  Co., 

New  York. 

Meyer,  R.E.  1980  Exponential  asymptotics.  SIAM  Rev.  22,  213-223. 

Mied,  R.P.  Sl  J.P.  Dugan  1975  Internal  wave  reflection  by  a  velocity  shear  and  density  ano¬ 
maly.  J.  Phys.  Oceanog.  5,  279-287. 

Miles,  J.W.  1961  On  the  stability  of  heterogeneous  shear  flows.  J.  Fluid  Mech.  10,  496-508. 

Miles,  J.W.  1963  On  the  stability  of  heterogeneous  shear  flows.  Part  2.  J.  Fluid  Mech.  16, 
209-227. 

Munk,  W.H.  1981  Internal  waves  and  small-scale  processes.  In  Evolution  of  Physical  Oceanog¬ 
raphy,  8.A.  Warren  and  C.  Wunsch  eds.,  MIT  Press,  264-291. 

Mysak,  L.A.  1976  A  kinetic  theory  for  internal  waves  in  a  randomly  stratified  fluid.  Dyn. 

Atoms.  Oceans  l  91-98. 

Olver,  F.W.J.  1974  Asymptotics  and  Special  Functions.  Academic  Press,  New  York. 

Phillips,  O.M.  1966  The  Dynamics  of  the  Upper  Ocean.  Cambridge  University  Press,  London 
and  New  York. 

Pinkel,  R.  1981  On  the  use  of  Doppler  sonar  for  internal  wave  measurements.  Deep-Sea  Res. 

28,  269-289. 

Sanford,  T.B.  1975  Observations  of  the  vertical  structure  of  internal  waves.  J.  Geophys.  Res. 

80,  3861-3871. 


Snow,  C.S.  1942  The  Hypergeometric  and  Legendre  Functions  with  Applications  to  Integral 


47 


Equations  of  Potential  Theory.  National  Bureau  of  Standards. 

Thorpe,  S.A.  1969  Neutral  eigensoluticins  of  the  stability  equation  for  stratified  shear  flow.  J. 
Fluid  Mech.  36,  673  -683. 

Thorpe,  S.A.  '  ”78  On  internal  gravity  waves  in  an  accelerating  shear  flow.  J.  Fluid  Mech.  88, 
623-639. 

Thorpe,  S.A.  1981  An  experimental  study  of  critical  layers.  J.  Fluid  Mech.  103,  321-344. 

van  Duin,  C.A.  A  H.  Kelder  1982  Reflection  properties  of  internal  gravity  waves  incident 
upon  a  hyperbolic  tangent  shear  layer.  J.  Fluid  Mech.  120,  50S-S21. 


APPENDIX  A. 


THE  REFLECTED  AND  TRANSMITTED  WAVE  AMPLITUDES 
IN  TERMS  OF  HEUN  FUNCTIONS. 


In  chapter  four,  it  is  shown  that  the  the  solution  corresponding  to  the  reflected  wave  as 
2  —  -oo  is  expressed  in  terms  of  the  Heun's  function  as 

h»i  -  F(ci,b-,a,0,y,8;U\)  .  (Al) 

This  is  the  solution  in  (4.14)  about  U i  «  0  with  zero-exponent.  The  solution  with  exponent 
— 2i>!  corresponds  to  the  incident  wave  as  z  —  and  can  be  written  in  terms  of  Heun's 
function  using  the  transformation  w2  —  U\  2<M,W|.  From  Snow  (1942), 

w2  ”  Ui  lli{F(.Ci,b2,l+or-y,l+l3—y,8.,U\)  ,  (A2a) 


where 

b2  —  b  —  (1 — y){S+C|(l-f-5+/9 — y— 8)}  .  (A2b) 

The  solution  belonging  to  the  zero-exponent  at  U\  •  1  in  (4.14)  corresponds  to  the 
transmitted  wave  as  z  —  -H».  Transforming  U i  —  \—U\  leads  to  the  solution  (Snow,  page 
123) 

wj  ™  F(l— C|,— b— c?/3;a,/8,l+a+/3 — y— 8,8;1— U\)  .  (A3) 


To  obtain  the  amplitudes  of  the  reflected  and  transmitted  waves,  the  solutions  w,  and 
w2  must  be  related  to  the  solution  w3.  For  the  Hypergeometric  functions,  these  connection  for¬ 
mulae  can  be  expressed  in  terms  of  gamma-functions,  but  for  Heun's  functions,  it  is  necessary 
to  match  the  solutions  and  their  derivatives  at  points  between  the  singularities.  The  matching 
point  must  be  positioned  at  a  distance  less  than  (cj — 1)  away  from  U\  —  1  to  be  within  the 
radius  of  convergence  of  w3.  If  we  denote  this  position  by  x,  it  follows  that 


and 
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y»'2w2  +  w2»>'2 

H»']  W2  +  Wj  Wj 
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(A4) 
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w2w\  -  w'jWi 
(Vjw'l  +  w'jW| 


(A5) 
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where  a  prime  indicates  differentiation  with  respect  to  U\  and  the  subscript  x  implies  that  all 
functions  are  evaluated  at  U\  -  x. 

The  numerator  in  (AS)  is  the  negative  of  the  Wronskian  for  the  two  solutions  W|  and 
w2,  which  we  evaluate  as 

"i 

-T  m 

Wjw'j  -  w',w2  -  °  (A6) 

where  p  is  the  coefficient  of  wU{  term  in  Heun’s  equation  (4.8).  Substituting  for  p  leaves 

_ 1 — y _  _ 1 _ 

ti 

The  same  procedure  can  be  applied  to  the  case  involving  a  critical  layer  within  the 
shear  region,  although  a  branch  cut  must  be  chosen.  A  lengthy  calculation,  which  will  not  be 
given  here,  leads  to  the  results  presented  in  figure  A1  and  in  Munk  (1980)  for  c— 0.  (Van 
Duin  &  Kelder  (1982)  give  a  simpler  technique  for  determining  the  reflection  and  transmission 
coefficients.)  Waves  are  overreflected  when  |/t|2  >  1,  which  requires  Ri  <  .13. 

For  values  in  the  unstable  region  of  the  a—J  plot  of  figure  A1  (Drazin,  19S8),  an 
infinite,  reflected,  intrinsic  energy  density  is  obtained.  This  is  the  case  of  resonant 
overreflection.  According  to  linear  theory,  the  shear  layer  radiates  waves  without  being  forced. 
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Figure  Al.  Reflection  from  a  critical  layer.  Plotted  are  the  values  of  the  fractional  intrinsic 
energy  density  reflected  and  transmitted  through  a  mean  shear  flow,  f/-£/0  tanh 0z  as  a  func¬ 
tion  of  minimum  Richardson  number  The  plot  is  drawn  for  /-2a2,  where 

a— kfi  is  the  dimensionless  horizontal  wavenumber.  This  corresponds  to  a  wave  packet  travel¬ 
ing  at  an  inclination  of  45°  as  ?— ±<».  (Ri  -  a1  corresponds  to  the  limiting  case  of  vertical 
group  velocity  as  r— ±».) 


APPENDIX  B. 


CHANGES  IN  PACKET  AMPLITUDE  DUE  TO  DISPERSION  AND  GEOMETRICAL  SPREADING. 

Using  ray  theory,  we  derive  equations  that  describe  the  dispersion  and  geometrical 
spreading  of  a  localized  wave  packet  in  a  uniform  medium  of  n  dimensions.  The  derivation 
generalizes  the  one  presented  in  Section  C  of  Chapter  5  and  leads  to  the  stationary-phase  result 

A  -  r*  (Bl) 


for  large  times. 

In  a  uniform  medium,  the  ray  equation  for  the  action  density  A  takes  the  form 


dA  _  ,  a2n  &kj 

dt  bkjbkj  bx, 


(B2) 


(Hayes,  1970),  where  the  subscripted  variables  are  components  of  K  and  x,  repeated  indices  are 
summed  over  ail  values,  and  the  dispersion  relation  is 


«-n(K(x,r))  . 


(B3) 


(For  convenience  in  notation  we  do  not  use  subscripts  for  partial  derivatives  in  this  appendix.) 
To  obtain  bk,/bxj,  we  write 

K  -  K(x  -  C tt)  ,  (B4) 

appropriate  to  a  uniform  medium,  and  define 


a  —  x  —  Crt 

Differentiating  (B4)  with  respect  to  x  leaves 


bki  Ski  ba , 
bxj  bap  bxj 


,  bk,  J  bk,  Va  a*„] 

baj  ba,  bkpbkm  bxj  j 

where  bkjbaj  is  the  value  of  bk,/bxf  at  t  -  0.  We  rewrite  (B6)  as 


(B5)  ' 


(B6) 
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which  in  the  long-time  limit  becomes 


Bkm_  i*L  Q2n  1  d*/ 

Bxj  ( Bap  Bkp  Bkm  j  t  Bctj 

Multiplying  (B8)  by  BaJBk,  and  using 

Baq  Bk, 

Bk,  Bot„ 


(B8) 


(B9) 


leaves 


Mm  B2  n  , 

Bxj  BkpBkm  “  t  ‘ 


(BIO) 


We  use  the  relation  8y>  -  n,  where  n  is  the  number  of  spatial  dimensions,  and  set  q  -  j  in 
(BIO)  to  obtain 


Bkm  ^  rt_ 

Bxj  Bkj  Bkm  t 

substituting  (Bll)  into  the  action  equation  (B2)  gives 

dA  .  n 

~dF~~A  7 


(BID 


(B12) 


or 


At"  —  constant  , 


(B13) 


which  matches  the  stationary-phase  result  (BD. 

We  note  the  following: 

(i)  Equation  (B7)  applies  to  those  components  of  K  that  lie  in  a  direction  in  which  the 
medium  is  homogeneous,  even  if  the  medium  is  not  homogeneous  in  ail  directions. 

(ii)  Not  all  initial  conditions  lead  to  the  result  given  by  (B13).  If,  in  (B7), 


-Jin  ii.  <0 

BkpBkm  Bap 


(B14) 


then  the  term  in  parentheses  in  (B7)  can  become  zero,  giving  an  infinite  value  for 
BkJBxj  and  for  the  action  density.  Physically,  the  Fourier  components  with  slow  group 


velocities  precede  those  with  faster  group  velocities  at  the  initial  time,  and  as  the  fast 
waves  overtake  the  slow  ones,  a  caustic  develops. 
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(iii)  We  list  here  for  use  in  Section  C  of  Chapter  5  the  components  of  d2Cl/(dk,dkj).  The 
dispersion  relation  is 


n-rrrn 

1  kf  +  ki  +  m2  J 

(B15) 

-HI 

a 

(B16) 

where  K  -  (ku  k2,  m)  and 

a  S  (ki  +  k})N2  +  m2/2 

(B17) 

b  =  (kf  +  kj  +  m2)  . 

(B18) 

Then  for  /  -  1,  2, 


dHt 

dk,2 


m2(N2  -  1  -  +  N2a~') 


a2n 

dM*2 


-m2(N2  -  f2)(ab)~V2  k\k2(la  +  MV2) 


(B19) 

(B20) 


d2U 

Cx,dm 


nk,(N2  -  f2)(a~v‘b~i/2)  2  -  m2(3b~'  +  /V) 


(B21) 


t  > 


with  d2fl/dm2  given  by  (5.10). 


APPENDIX  C 


WAKES  AND  MEAN  FLOWS  OF  LOW-FREQUENCY  INTERNAL- WAVE  PACKETS* 


A.  Introduction 

The  role  of  wave  transience  in  affecting  mean  flows  and  mean  profiles  of  temperature 
and  tracers  is  well-established,  a  prime  example  being  the  sudden  stratospheric  warming,  where 
growing  planetary  waves  modify  the  temperature  distribution  by  as  much  as  40°C  and 
decelerate  upper  level  winds. 

Transience  is  a  marked  feature  of  energetic,  near-inertial  oscillations  in  the  ocean. 
Although  it  is  unlikely  that  they  will  cause  ocean  temperatures  to  increase  by  40°C,  these  oscil¬ 
lations  could  be  very  effective  at  inducing  mean  flows.  If,  as  for  high-frequency,  plane  internal 
wave  packets, 

u-uL  -  0(4)  ,  (Cl) 

c 

where  u  and  u£  are  the  Eulerian  and  Lagrangian  horizontal  mean  velocities,  respectively,  £  is 
the  intrinsic  energy  density  (see  equation  (C16)  and  below),  and  c'  is  the  horizontal,  intrinsic 
phase  velocity,  then  for  fixed  horizontal  wavelength,  the  generation  of  mean  velocities  would 
be  heavily  dominated  by  low-frequency  waves:  in  addition  to  the  phase  velocity  appearing  in 
the  denominator  in  (Cl),  energy  spectra  typically  fail  off  as  the  inverse  square  of  the  frequency 
and  even  faster  near  the  inertial  and  tidal  bands.  A  rough  calculation  using  (Cl)  and  the  obser¬ 
vation  that  low-frequency  waves  are  marginally  stable  (i.e.  have  a  minimum  Richardson 
number  approaching  unity),  gives  mean-flow  magnitudes  of  nearly  half  the  orbital  velocity  of 
the  waves. 

When  Coriolis  terms  are  important,  however,  mean  flows  can  evolve  differently.  Here 
mean-flow  equations  are  derived  for  low-frequency  internal  wave  packets.  Of  the  cases  con¬ 
sidered,  only  when  rotation  is  ignored  does  u  or  uL  equal  E/c  identically,  and  only  when  hor¬ 
izontal  dependence  is  ignored  does  an  inertial  wave  appear  in  the  wake  of  the  packet,  as 
predicted  by  Hasselmann  (1970)  for  a  wave  field  that  is  statistically  homogeneous  in  the  hor¬ 
izontal.  These  calculations  supplement  the  related  work  of  Andrews  (1980),  Bretherton  (1969), 
Grimshaw  (197S),  and  Hasselmann  (1970). 

* 5tefer*ncei  ire  included  at  the  end  of  this  Appendix. 
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It  is  assumed  that  the  VaisaJa  frequency  N  is  constant  and  that  there  is  no  background 
current.  For  the  upper  thermocline,  where  low-frequency  packets  are  most  energetic  and 
identifiable,  these  are  poor  assumptions,  but  they  help  to  isolate  the  dependence  of  mean-flow 
evolution  on  wave  frequency.  The  calculations  can  probably  be  extended  to  include  variations 
in  the  medium  through  a  WKB  analysis  without  altering  the  conclusions  obtained  below.  We 
also  ignore  ocean  boundaries,  which  is  reasonable  for  the  scales  to  be  considered:  vertical 
packet  dimensions  of  a  few  hundred  meters,  implying  horizontal  wavelengths  of  roughly  S-20 
kilometers,  depending  on  frequency. 

The  equations  to  be  solved  are 

n,  +  n-Vu  +  2flxu  +  — Vp  +  g  —  0  (C2a) 

Po 

V  u  -  0  .  (C2b) 

The  fluid  is  nondissipative,  Boussinesq,  and  incompressible,  211  -  (0,0,/),  where  /  is  the  iner¬ 
tial  frequency,  assumed  constant,  and  a  —  (u,v,  w). 

B.  Mean  flows 

As  waves  propagate  through  a  fluid,  they  induce  stresses  that  can  accelerate  mean 
motions.  There  are  many  examples  with  physical  interpretations  in  the  literature.  Mathemati¬ 
cally,  a  mean  flow  arises  because  a  sinusoid  is  generally  not  an  exact  solution  to  the  equations 
of  motion.  If  sinusoidal  motion  is  assumed,  the  nonlinear  terms  generate  harmonics  and  mean 
quantities.  Even  in  the  special  cases  where  sinusoidal  motion  is  an  exact  solution  (e.g.  plane 
internal-inertial  gravity  waves  in  a  uniform,  nondissipative,  Boussinesq  fluid)  strictly  sinusoidal 
motion  may  not  be  consistent  with  the  condition  of  no  motion  initially.  For  instance,  a  wave 
packet  with  a  slowly  growing  leading  edge  induces  mean  as  well  as  orbital  motion.  Where  the 
amplitude  of  the  wave  packet  is  constant,  the  horizontal  velocity  is  ii  -  *sin(Kx-o»r)  +  u, 
where  u  is  an  0( t2)  constant.  This  is  still  an  exact  solution  under  the  assumptions  mentioned 
above  provided  that  u  is  the  Doppler-shifted  frequency,  related  to  the  intrinsic  frequency  <»’  by 
oi  —  «'  +  Uk  for  k> 0. 

When  rotation  (and  dissipation)  can  be  ignored,  net  accelerations  produced  by  the 
leading  edge  of  a  wave  packet  are  cancelled  by  decelerations  in  the  trailing  edge.  The  mean 
flow  returns  to  zero  after  the  packet  passes,  although  a  net  displacement  of  the  fluid  occurs  if 
the  Lagrangian-mean  flow  is  nonzero. 


When  rotation  is  important,  the  trailing  edge  may  not  decelerate  the  fluid  back  to  zero 
velocity  because  the  mean  motion,  now  affected  by  the  Coriolis  force,  changes  direction  during 
the  passage  of  the  packet.  The  same  forces  that  set  up  the  motion,  therefore,  may  no  longer 
stop  it.  Instead,  the  system  is  left  ringing  at  inertial  frequency  after  the  packet  propagates  away. 
This  is  illustrated  in  the  following  example. 

Solutions  to  (C2)  are  sought  of  the  form 

(u,v,w,p,9)  -  (u',v',w',p',9')  +  (u,v,w,p,9)  .  (C3) 

The  buoyancy  parameter  is  9  *■  gp'/po  and  a  factor  of  1/po  is  included  in  the  pressure  variable 
for  simplicity  in  notation.  The  primed  quantities  are  0(c)  and  behave  like  sin(K-x  -  o>r), 
where  K  *■  (k,Q,m).  The  barred  quantities  are  obtained  by  an  Eulerian  average  over  phase  and 
are  0(e2). 

The  0(e)  equations. 


u',  -  fv'  +  p'x  -  0 

(C4) 

V',  +  /«'  +  P'y  «  0 

(C5) 

w',  +  V  +  p\  -  0 

(C6) 

9\  -  tfV  -  0 

(C7) 

u\  +  v'y  +  w',  -  0  , 

(C8) 

yield  the  linear-theory  results:  the  dispersion  relationship 

(C9) 

and  the  relationships 

_ _ _  _  f  /«'  V*  i(N2—u>2)k  N*k) 

/m, 

'  (CIO) 

(Andrews,  1980). 

The  0(e*)  equations  describe  the  evolution  of  means  and  second  harmonics.  The 
later  component  is  eliminated  by  averaging,  leaving 


puT  -  Jv  +  ppx  -  -ptiPbx  -piu^z 


(Cll) 


fiwT  +  9  +  npz  -  -liOTw^x  -  fi(w'2)z 


(C13) 


H0T  ~  N2w  -  -M( v'O')  r  (C14) 

fl(Ux  +  Vy  +  WZ)  —  0  .  (CIS) 

These  equations  are  written  in  terms  of  the  stretched  variables  X  -  pix,  r-/*'  'V  «  1) 
which  vary  by  an  0(1)  amount  over  the  scales  of  the  packet.  The  boundary  condition  specifies 
that  all  quantities  vanish  as  Z  —  — «>  (envisioning  a  packet  with  downward  group  velocity). 
The  0(«)  equations  have  been  used  to  eliminate  products  of  quantities  on  the  right-hand-side 
of  (C11MC1S)  which  are  out  of  phase.  The  mass  source  in  the  Lagrangian-mean  continuity 
equation  (CIS)  does  not  appear  for  this  problem  until  0(e4). 

Case  (1)  3/3*  -  3/3  Y  -  0,  frp  -  0(1) 

By  ignoring  variations  in  X  and  Y,  the  continuity  equation  and  boundary  condition 
force  w  to  be  zero  everywhere.  The  buoyancy  equation  then  gives  3-0,  leaving  the  pressure 
gradient  term  to  balance  the  Reynolds  stresses  in  the  vertical  momentum  equation. 

A  high  frequency  approximation,  o'/f  »  1,  reduces  the  equation  for  the  x- 
component  of  the  pseudomomentum,  pu  to 

P\t  “  —  (ulXDz  (C16) 

where  to  Of*2),  p\  —  Etc '  and  E  —  (u^+  v^+  w^  +  3^/A*)/2,  and  c'  —  o' Ik.  The  horizon¬ 
tal  momentum  equations  then  become 


puT-fv-nP\T  (C17) 

M vt  +  fu  -  0  .  (C18) 

At  this  stage  we  assume  that  the  Coriolis  terms  and  the  acceleration  terms  are  the 
same  order,  i.e.  that  frp  —  0(1),  where  rp  is  the  time  scale  of  the  packet.  (Since  rp  »  I/o' 
is  necessary  for  scale  separation,  this  assumption  includes  the  high  frequency  approximation 
leading  to  (06)).  The  case  of  strong  rotation,  frp  »  1  is  considered  in  the  next  section. 

Defining  4  —  u  +  (v,  (07)  and  (08)  combine  to  become 

p$T  —  ifa  —  pP\t  »  (09) 

which  has  the  solution 

5  -  j  -  ifei(}l»)T  f  e-{//»)r£(r)dr  .  (C20) 
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The  second  term  arises  when  /  is  nonzero  and  represents  an  inertial  wave  left  in  the  wake  of 
the  packet.  Hasselmann  (1970)  obtains  a  similar  solution  for  a  random  field  of  surface  waves 
or  internal  waves  in  shallow  water.  Related  behavior  illustrating  that  Rossby-wave  packets 
leave  lower-frequency  Rossby  waves  in  their  wake  is  shown  in  Rhines  and  Holland  (1979). 

It  therefore  appears  that  transient  waves  need  not  induce  transient  mean  flows,  even  in 
a  nondissipative  fluid.  This  conclusion  follows  more  directly  when  the  equations  are  rewritten 
in  Lagrangian-mean  form: 


-  fvL  +  Px  “  ~p(pDx  -  p(p^Y~  M(pnPz 

"PP\T 

(C21) 

litf  +  fuL  +pLr‘m  0 

(C22) 

+  9L  +  pL  Z-n(p£t)x  -  r)  Y  -  M^PnPz 

”  ~PPiT 

(C23) 

o 

1 

r* 

1 

1^ 

4. 

(C24) 

piufl  +  vf  +  w£)  —  0 

(C25) 

where  f(x,t)  *»  is  the  particle  displacement,  and  the  psuedomomemtum  vector  p  has 

components  Cpi,0,p3).  The  Lagrangian  mean,  defined  as 

g(x,t)L  -  g(x+f,t)  (C26) 

for  any  quantity  g ,  is  an  ensemble  average  over  phase.  It  is  then  possible  to  define  f  such  that 
{ (x,t)  —  0  identically  (Andrews  St  McIntyre  (1978)).  Dependence  on  X  and  Y  in  (C21)- 
(C25)  has  been  included  for  later  reference.  Ignoring  these  dependences  for  now,  the  X- 
momentum  equation  can  be  rewritten  as 

pXuL  -pOt-/vl  -0,  (C27) 

which  in  turn  can  be  rewritten  as 

DL  (uL  -  px  -  fy)  -  0 

(Andrews  St  McIntyre  (1978)).  DL  signifies  the  rate  of  change  moving  with  the  Lagrangian- 
mean  flow: 
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DL-&+iLix  >cm 

If  the  Coriolis  force  leads  to  a  net  change  in  y  during  the  passage  of  the  packet,  uL  must 
remain  nonzero  after  the  packet  propagates  away. 

Case  (2)  B/BX  -  B/BY-0  ft,  »  1 

The  solution  given  in  (C20)  is  derived  assuming  that  the  Coriolis  terms  are  small  but 
nonzero.  This  requires 

(C2„ 

We  identify  r„  a  T/ft  with  the  time  scale  of  the  packet.  Let  tp  -  2 nn/oi'  with  n  measuring 
the  number  of  wave  periods  required  for  the  wave  packet  to  reach  maximum  amplitude  at  a 
fixed  location.  Then  frf  —  0(1)  translates  to  «'  —  Oilirnf).  If  n  «■  4,  say,  to'  must  be  in 
the  range  of  25/  for  the  Coriolis  and  acceleration  terms  to  be  of  the  same  order.  Since  25/ 
corresponds  to  a  wave  period  of  less  than  one  hour  at  midlatitudes,  this  division  lies  above  the 
frequency  of  the  semidiurnal  tide  and  tidal  harmonic  at  six-hour  period.  The  assumption 
ftp  —  0(1)  therefore  excludes  what  is  by  far  the  most  energetic  part  of  the  oceanic  internal 
wave  spectrum.  The  amount  that  is  excluded  depends  of  course  on  n,  but  the  point  is  that  it  is 
hrnf/u  and  not  //«'  that  measures  the  ratio  of  the  Coriolis  to  acceleration  terms  in  the  mean 
equations. 

The  approximate  balance  to  lowest  order  is  now  between  the  Coriolis  terms  and  the 
radiation  stresses,  so  that 


~*J  P\T 

(C30) 

U }  . 

JfPlTT- 

(C31) 

Because  the  mean  motion  is  reduced  to  an  Oipi)  quantity,  the  Stokes  drift  appears  at  the  same 
order 

vs  -  (Wv’)z  ■  (C32) 

More  interestingly,  vL  and  uL  return  to  zero  after  the  packet  propagates  away.  If  the  mean 
velocities  are  expended  in  powers  of  p,  inertial  motion  in  the  wake  still  does  not  appear,  for  at 
all  higher  orders  vL  and  uL  can  be  related  to  higher-order  derivatives  of  px  with  respect  to  T. 


To  see  what  has  happened  to  the  inertiaJ  motion,  we  can  derive  an  exact  solution  for 
<frL  -  uL  +  tvL  analogous  to  (20)  for  <t>  except  that  when  written  in  terms  of  Lagrangian-mean 
variables,  the  solution  is  valid  for  all  frequencies:  no  high-frequency  approximation  is  neces¬ 
sary  to  replace  the  right  hand  sides  of  (C21)  and  (C23)  with  p\T  and  piT  respectively.  (Pseu¬ 
domomentum  is  not  exactly  conserved  for  localized  wave  packets  since  mean  quantities  depend 
on  position;  however,  corrections  do  not  appear  in  the  present  context  until  0(e4).)  We  can 
thus  write 

!>L  -  j e‘^,ltHT-nfipirdr  .  (C33) 

Since  by  assumption  frt  »  1,  the  phase  of  the  exponential  function  in  the  integrand  of 
(C33)  varies  rapidly.  An  asymptotic  expansion  can  then  be  found  through  an  integration  by 
parts.  Integrating  N  times  leaves 

lL  -  -^Pxt  +  •  ■  •  +  f  e^^j—dr  .  (C34) 

The  last  integral  represents  the  inertial  motion  left  in  the  wake  of  the  packet.  It  does  not 
appear  in  the  asymptotic  expansion  because  it  is  exponentially  small  as  n  —  0.  (Using  a 
stationary-phase  calculation  and  noticing  that  the  phase  function  in  the  integral  has  no  station¬ 
ary  points,  the  integral  can  be  shown  to  be  exponentially  small  as  fi  vanishes  (Lighthill,  1978).) 


Case  (3)  b/~>X,  b/bY  ~  0(1);  /r,  »  1 

Low-frequency,  internal  wave  packets  have  horizontal  scales  that  are  much  greater 
than  their  vertical  scales;  however,  unlike  Hasselmann  (1970),  we  consider  horizontal  varia¬ 
tions  because  the  ratio  of  the  Reynolds  stress  terms  in  the  x-momentum  equation,  for  example, 
is 


jn  A Z 

(«V)2  k  A*  ’ 


(C35) 


where  A Z  and  &X  are  the  vertical  and  horizontal  scales  of  the  packet.  The  ratio  thus  scales  as 
the  number  of  horizontal  wavelengths  to  the  number  of  vertical  wavelengths,  not  as  the  ratio  of 
absolute  packet  dimensions  or  wavenumbers.  A  similar  argument  applies  to  the  radiation 
stresses  in  the  Lagrangian-mean  equations. 


Allowing  for  the  packet  to  produce  mean  horizontal  pressure  gradients  and  allowing  for 
all  the  terms  in  the  continuity  equation  leaves  the  system  of  equations  given  in  (C2IHC25). 
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We  write 


yi  _  y(0>  +  +  pt2  W<2)  +  •  •  ■  ,  (C36) 

and  similarly  for  ail  other  mean  quantities,  and  attempt  an  expansion  which  balances  terms 

geographically  and  hydrostatically  at  lowest  order: 

-/v<®  +  pi®  -  0  (C37) 

/tt(0)  +  pj®  -  0  (C38) 

*«»  +  pp  »  0  (C39) 

wi®  +  vi®  -  0  .  (C40) 


A  stream  function,  can  be  defined  at  this  order,  related  to  p(0>  by 

-  VAJpl®  ,  (C41) 

where  Vj|  -  a2/a*2  +  32/d  Y2,  u(®  -  (?>,  and  v(®  -  * #». 

The  evolution  of  the  0(1)  quantities  is  obtained  from  the  0(p)  equations: 


«i®  - /v(,)  +  piu  -Pit 

(C42) 

vi®  +  /«(,)  +  piu  »  0 

(C43) 

*4® +  «<■>  + pi” -p3r 

(C44) 

gjP>  _  jyV”  _  0 

(C45) 

«i”  +  vf”  +  **i”  —  0 

(C46) 

These  can  be  manipulated,  using  (C41),  to  obtain  Poisson’s  equation  for  </»(® 

^bl>(0)  +  ^^--Pu  ■  (C47) 

The  formal  solution  to  (C47)  can  be  written  in  terms  of  a  Green’s  function  or  derived  using  a 
three-dimensional  Fourier  transform. 


Suppose  we  choose  a  Gaussian  packet  envelope: 

Pi  ~  exphaj'tf-C.in2  -  a}Y2  -  a/(Z-C„ D2) 


(C48) 
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where  Cf  —  (C,i,0,Cfj)  is  the  group  velocity.  Symmetry  considerations  lead  to 

-  pP  -  -  P?')  r-o  -  0  ,  (C49) 

so  v4  is  zero  along  the  T-centerline  of  the  packet,  where  the  wave  amplitude  reaches  a  max¬ 
imum. 

In  the  far-field,  <£(o)  and  pM  decay  as  l/r2  multiplied  by  some  angular  dependence. 
(An  electromagnetic  analogy  is  the  potential  due  to  a  dipole  source.)  Consequently,  no  inertial 
wake  results.  Viewed  another  way,  (47)  becomes  Laplace’s  equation  as  T  — -  «>,  X  fixed, 
which,  because  of  the  boundary  conditions,  cannot  support  motion. 

As  in  Case  (2)  the  inertial  motion  does  not  appear  at  higher  order.  The  Oin2)  equa¬ 
tions  give 

-  Pvcr  +  jjiPiZT  .  (C50) 

which  also  reduces  to  Laplace’s  equation  after  the  packet  propagates  away.  Similarly,  the  Oin}) 
equations  lead  to 

l+$  •  (C51) 

Since  the  structure  of  the  equation  at  all  higher  orders  is  the  same  as  (C51),  the  expansion 
yields  only  decaying  motion  in  the  packet  wake. 

Again  as  in  Case  (2),  the  inertial  motion  may  result  from  an  expression  which  is 
exponentially  small  as  n  —  0.  There  is  at  least  one  other  possibility  though:  the  inertial  wake 
may  be  limited  to  packets  of  infinite  horizontal  extent,  i.e.  the  limit  of  B/BX,  d/d  Y  approaching 
zero  may  be  a  singular  one.  The  idea  is  that  although  the  horizontal  pressure  gradients  build 
slowly  in  a  long  but  finite  packet,  they  can  eventually  reach  an  0(1)  magnitude.  The  dynamics 
change  when  horizontal  pressure  gradients  are  important  and  inertial  oscillations  may  disappear. 

An  exact  solution  to  (C21MC25)  would  reveal  whether  the  limit  B/BX,  d/d  T  —  0  is 
singular.  A  suggestion  made  by  Glenn  Flieri  (personal  communication)  is  to  try  a  one-cycle 
sinusoidal  packet 

Pi  —  sink  ,  X  -  (a-X-a  T)  ,  (C52) 

for  ~ it  <  X  <  ir  and  p\  -  0  otherwise.  This  allows  solutions  to  be  obtained  by  assuming 
behavior  of  the  form  /4cosX  +  fisinx  for  each  of  the  mean  variables.  To  meet  the  initial 
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condition,  solutions  to  the  homogeneous  equations  are  also  needed.  The  effect  of  discontinui¬ 
ties  in  77  at  X  -  —it, it  requires  assessment. 

Another  possibility  is  to  ignore  Z  dependence.  The  equations  then  reduce  to 


ntf- fvL  +  npjt  -  np\  T  (C53) 

nvj>  +  fuL  +  pp{>  -  0  (C54) 

m(«>  +  v£)  -  0  ,  (C55) 

which  can  be  combined  to  give  the  "wake-less”  equation 

V^-p,k;  (C56) 


however,  if  X  and  Y  variations  are  ignored  from  the  outset,  so  the  forcing,  p\T,  builds  in  time 
everywhere  by  the  same  amount,  then  the  governing  equation  becomes  the  "wake-full"  equa¬ 
tion  (C19).  This  model,  although  somewhat  contrived,  is  an  example  illustrating  the  singular 
limit  associated  with  the  neglect  of  horizontal  dependences. 


C.  Non  uniform  metis  and  changes  in  packet  amplitude 

Grimshaw  (197S)  derives  a  general  system  of  mean  flow  equations  that  allows  for  vari¬ 
ations  in  the  Vaisala  frequency  and  for  background  shears.  Even  in  a  uniform  medium,  how¬ 
ever,  the  wave  induced  mean  flow  will  be  affected  by  changes  in  packet  amplitude  resulting 
from  dispersion  and  geometrical  spreading,  especially  when  the  propagation  is  three- 
dimensional  and  especially  when  the  waves  are  highly  transient  or  localized  and  therefore 
involve  a  broad  (or  at  least  non-narrow)  range  of  frequencies  and  wavenumbers.  A  stationary- 
phase  calculation  describes  the  attenuation  of  energy  in  a  uniform  medium,  but  this  result  can 
also  be  obtained  from  the  ray  equations.  The  equation  for  action  density  along  the  ray  (in  one 
dimension  for  now)  is 


J£) ,  (C57) 

where  u(x,t)  -  (l(k,x,t)  is  the  dispersion  relation  and  d/dt  -  3/3/  +  C,  B/dX.  Since 
<*/<*- 0 


kx  -  [k(x-Ctt))x 


o 

i+n**(xT)  ,_of 


(C58) 


Substituting  (C58)  into  (CS7)  and  using  —  0,  gives  the  stationary  phase  result  that  the 
action  density  decays  as  1/r  as  t  — >  «>. 

It  may  seem  surprising  that  although  dispersion  results  from  a  spread  in  frequencies  or 
wavenumbers,  bandwidth  information  does  not  enter  directly  into  the  equation  for  action  den¬ 
sity,  as  it  does  in  the  equation  for  action  spectral  density  n(k,x,t ): 
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bn  .  bn  dk  .  bn  dx 
bt  bk  dt  bx  dt~ 


(C59) 


Here  spectral  content  appears  through  the  term  bn/bk\  recall,  however,  that  the  bandwidth 
does  not  enter  in  the  stationary  phase  calculation  either,  except  at  higher  order  when  developed 
as  an  asymptotic  expansion  (Fdsen  and  Marcuvitz,  1973).  Writing 

oo 

f(x,t)  -  f  g(k)e‘*lk)'dk  ,  (C60) 


where  <it(k)  —  kx/t  —  <n(k),  we  transform  variables,  using 

J2  -  /M*c)  -  ifr(Ar)) 

(Note  that  using  L’Hopital's  rule, 

G(0)  -  gikJ^J 

provided  +”(ko)  is  nonzero.)  Equation  (C60)  then  becomes 

f  -  (G( 0)  +  (7'(0)s  +  C"(0)y  +  •  •  •  )e~l2,ds  . 


(C61) 

(C62) 


(C63) 


(C64) 


The  stationary  point  is  k^,  which  satisfies  t/»’Oco)  -  0,  and  a  prime  in  (C63)  denotes 
differentiation  with  respect  to  s.  The  bandwidth  information  in  G'(0),  G"( 0)  •  •  •  appears  only 
as  a  correction  to  the  iowest-order  solution. 

The  advantage  of  the  ray  equation  (CS7)  over  the  method  of  stationary  phase  is  that  it 
is  easily  extended  to  non  uniform  media.  Complications  arise  because  in  general  kx  is  not 
known  analytically  and  because  in  three  dimensions,  nine  equations  (six  of  which  are  indepen¬ 
dent)  are  required: 


-v»uu 


(C65) 
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(Hayes,  1970).  Still  the  system  of  equations  can  be  integrated  numerically.  As  an  approxima¬ 
tion,  it  is  probably  justified  to  ignore  horizontal  variations  in  the  Vaisala  frequency  and  shear, 
so  that  the  horizontal  wavenumbers  are  conserved,  but  the  point  here  is  that  for  localized  pack¬ 
ets,  dispersion  and  geometrical  spreading  due  to  a  spread  in  horizontal  wavenumbers  should  not 
be  dismissed  as  unimportant.  Its  affect  on  packet  amplitudes  may  match  that  produced  by 
refraction. 


Conclusion 


The  character  of  the  mean  flow  induced  by  an  internal  wave  packet  is  sensitive  to  the 
frequency  of  the  waves  and  the  dimensions  of  the  packet. 


When  horizontal  variations  are  neglected,  high-frequency  waves  induce  a  mean  flow 
given  by  uL  —  E/c,  unless  rotation  is  included.  Then  uL  is  given  by  (C20)  and  an  inertial 
wave  oscillates  in  the  wake  of  the  packet.  As  the  wave  frequency  decreases  toward  f,  the  mean 
flow  adjusts  so  as  to  move  perpendicular  to  K  at  lowest  order,  moreover,  in  contrast  to  the 
high-frequency  case,  uL  vanishes  whenever  the  packet  amplitude  levels  off  (see  (C29)  and 
(C30)). 


When  horizontal  dependences  are  included,  thn  lowest  order,  mean  flow  equations  lead 
to  Poisson’s  equation  for  the  pressure  or  stream  function  regardless  of  wave  frequency. 
Motions  induced  by  low  frequency  packets  balance  hydrostatically  and  geostrophically.  A  ques¬ 
tion  suggested  by  these  calculations  is  whether  the  limit  B/BX,  d/6  Y  —  0  is  a  singular  one,  and 
thus  whether  it  is  realistic  to  expect  wakes  with  inertial  oscillations. 


Table  1  provides  a  summary  of  results. 
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